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Summary 


In  this  report  we  describe  the  application  of  a  modern,  time-accurate,  Euler/Navier- 
Stokes  analysis,  or  numerical  unsteady  aerodynamic  simulator,  to  the  prediction  of  unsteady 
subsonic  and  transonic  flows  through  a  two-dimensional  cascade.  In  particular,  we  consider 
unsteady  flows  excited  by  harmonic,  pure  torsional  and  pure  bending  vibrations  of  the  blades 
of  the  Tenth  Standard  Cascade.  One  purpose  of  this  study  is  to  validate  the  Euler/Navier 
Stokes  analysis  along  with  an  existing  linearized  inviscid  analysis  via  result  comparisons  for 
unsteady  flows  that  are  representative  of  those  associated  with  blade  flutter  in  modern  fans 
and  compressors.  We  have  also  applied  the  nonlinear,  time-accurate  analysis  to  determine 
the  relative  importance  of  nonlinear  and  viscous  effects  on  the  unsteady  response  of  a  cas¬ 
cade  to  prescribed  blade  motions.  It  is  anticipated  that  such  applications  will  guide  the 
further  development  of  time-accurate,  nonlinear  analyses  as  well  as  the  approximate,  but 
very  efficient,  unsteady  aerodynamic  analyses  that  are  used  in  turbomachinery  aeroelastic 
and  aeroacoustic  design  prediction  systems. 

The  results  of  the  present  study  reveal  a  close  agreement  between  inviscid  and  viscous 
blade  loadings  for  unsteady  subsonic  flows.  Also,  the  unsteady  surface  pressure  and  local 
work  per  cycle  responses  in  such  flows  are  essentially  linear,  and  they  are  predicted  quite 
accurately  using  a  linearized  inviscid  analysis.  For  unsteady  transonic  flows,  shocks  and 
their  motions  introduce  signiflcant  nonlinear  or  second  and  higher  harmonic  contributions  to 
the  local  unsteady  response.  Viscous  displacement  effects  tend  to  diminish  shock  strength 
and  the  impulsive  unsteady  loads  associated  with  shocks  and  their  motions.  The  energy 
transfer  between  the  fluid  and  structural  motions  in  both  subsonic  and  transonic  flows  is 
essentially  captured  by  the  flrst-harmonic  component  of  the  nonlinear  inviscid  and  viscid 
solutions,  but  in  transonic  flows  the  nonlinear,  flrst-harmonic  and  the  linearized  inviscid 
response  predictions  differ  significantly  in  the  vicinity  of  shocks  and  in  the  supersonic  regions 
upstream  of  the  shocks. 


1 


1.  Introduction 


A  time-accurate  Navier-Stokes  analysis,  or  numerical  simulator,  is  needed  for  understand¬ 
ing  the  relative  importance  of  nonlinear  and  viscous  effects  on  the  unsteady  flows  associated 
with  turbomachinery  blade  vibrations  and  blade-row  noise  generation.  Once  validated,  such 
an  analysis  could  also  serve  as  a  benchmark  for  evaluating  and  guiding  the  further  develop¬ 
ment  of  the  asymptotic  unsteady  aerodynamic  models  that  are  currently  used  in  turboma¬ 
chinery  aeroelastic  and  aeroacoustic  design  prediction  systems. 

Until  recently,  the  unsteady  aerodynamic  analyses  that  have  been  used  in  aeroelastic 
and  aeroacoustic  design  studies  were  based  on  classical  linearized,  inviscid,  flow  theory 
(see  [Whi87]  for  a  review),  which  essentially  applies  to  cascades  of  unloaded  flat-plate  blades 
that  operate  in  entirely  subsonic  or  entirely  supersonic  flow  environments.  Over  the  past 
decade,  more  general  inviscid,  unsteady,  aerodynamic  linearizations  have  been  developed 
(e.g.,  see  [Whi90,  Ver92,  HC93]).  These  account  for  the  effects  of  important  design  features, 
such  as  real  blade  geometry,  steady  blade  loading,  and  operation  at  transonic  Mach  numbers, 
on  the  unsteady  aerodynamic  response  of  blade  rows.  This  type  of  model  has  received  con¬ 
siderable  attention  in  recent  years,  and  numerical  methods  and  computer  codes  for  solving 
the  resulting  linearized  equation  sets  are  now  being  applied  in  turbomachinery  aeroelastic 
design  prediction  systems  [Smi90,  MM93]. 

Although  linearized  inviscid  analyses  meet  the  needs  of  designers  for  efficient  unsteady 
aerodynamic  response  predictions,  of  necessity,  they  ignore  potentially  important  physical 
and  geometrical  properties  of  the  flow,  including  the  effects  of  moderate  to  large  amplitude 
unsteady  excitation  and  the  effects  of  viscous-layer  displacement  and  separation.  Since 
the  early  1980’s,  a  number  of  time-accurate  Euler  and  Navier-Stokes  procedures  have  been 
developed  to  predict  unsteady  flows  (see  [Ver93]  for  a  review).  These  have  been  applied  to 
predict  flows  through  single  blade  rows  in  which  the  unsteadiness  is  caused  by  prescribed 
blade  vibrations  ([HR89,  He90,  Sid91])  or  by  prescribed  aerodynamic  disturbances  at  the 
inflow  or  outflow  boundaries  [Gil88],  and  flows  through  aerodynamically  coupled  arrays  in 
which  the  unsteadiness  is  caused  by  the  relative  motions  of  adjacent  blade  rows  ([Rai87, 
Rai89]). 

These  recent  and  important  advances  in  the  numerical  simulation  of  unsteady  turbo¬ 
machinery  flows  suggest  that  it  is  now  appropriate  to  carefully  validate  and,  if  necessary, 
extend  a  time-accurate,  Euler/Navier-Stokes  procedure  for  the  prediction  of  turbomachinery 
aeroelastic  and  aeroacoustic  response  phenomena.  Usually,  numerical  procedures  are  veri¬ 
fied  via  comparisons  with  experimental  data  (e.g.,  [MW92]),  but  because  of  the  numerous 
controlling  parameters  and  uncertainties  involved,  it  is  often  difficult  to  ascertain  causes  for 
the  differences  between  the  numerical  and  experimental  results.  In  the  present  effort,  we 
are  taking  an  alternative  approach;  one  in  which  solutions  based  upon  very  different  analyt¬ 
ical  procedures  are  compared,  both  to  validate  the  procedures  and  to  better  understand  the 
relevant  unsteady  flow  phenomena. 

A  validated  Navier-Stokes  analysis  for  turbomachinery  unsteady  flows  will  represent  an 
important  advance  in  unsteady  aerodynamic  theory.  This  numerical  simulator  (or  “wind 
tunnel”)  can  provide  engineers  with  useful  insights  on  the  nonlinear  and  viscous  effects 
associated  with  blade  vibration  and  discrete-tone  noise  generation.  It  can  also  provide  a 
test-bed  for  assessing  and  improving  the  asymptotic,  i.e.,  the  linearized  inviscid  and  high 


2 


Reynolds  number,  inviscid/viscid  interaction,  models  that  are  currently  being  developed  for 
use  in  turbomachinery  aeroelastic  and  aeroacoustic  design  prediction  systems. 

For  this  purpose,  in  a  previous  investigation,  an  existing  multi-blade- row  Navier-Stokes 
analysis,  i.e.,  the  R0T0R2  analysis  developed  by  Rai  [Rai87,  Rai89],  was  modified  and  ap¬ 
plied  to  predict  the  unsteady  flows  through  isolated,  two-dimensional  blade  rows.  In  this 
work  [DV94]  unsteady  subsonic  flows  excited  by  prescribed  aerodynamic  disturbances  at 
inlet  and  exit  that  carry  energy  towards  the  blade  row  (the  gust  response  problem)  were 
considered.  To  evaluate  this  single  blade  row  analysis,  called  NGUST,  numerical  simula¬ 
tions  were  performed  for  benchmark  inviscid  and  viscous  unsteady  flows,  and  the  predicted 
results  were  compared  with  those  based  on  linearized  inviscid  flow  theory.  For  small  am¬ 
plitude  unsteady  excitations,  the  unsteady  pressure  responses  predicted  with  the  nonlinear 
and  linearized  analyses  showed  very  good  agreement,  both  in  the  field  and  along  the  blade 
surfaces.  Based  upon  a  limited  range  of  parametric  studies,  it  was  also  observed  that  the  un¬ 
steady  responses  to  inlet  vortical  and  acoustic  excitations  were  linear  over  a  surprisingly  wide 
range  of  excitation  amplitudes,  but  acoustic  excitations  from  downstream  produce  responses 
with  significant  nonlinear  content. 

The  objectives  of  the  present  GUIde  program  are  to  further  develop  a  numerical  simula¬ 
tor  for  inviscid  and  viscous  unsteady  flows  through  two-dimensional  blade  rows,  and  to  gain 
insight  into  the  relative  importance  of  nonlinear  and  viscous  effects  on  unsteady  aerodynamic 
responses  to  prescribed  blade  motions  (the  flutter  problem).  In  particular,  time-accurate,  Eu¬ 
ler  and  Navier-Stokes  solutions,  based  on  the  NPHASE  analysis,  see  [SLH+94]  for  a  detailed 
description,  have  been  determined  for  subsonic  and  transonic  unsteady  cascade  flows  excited 
by  prescribed  pure  torsonal  and  pure  bending  vibrations  of  the  blading.  The  NPHASE  anal¬ 
ysis  which  is  based  upon  the  high-resolution,  wave-split,  implicit,  time  marching,  numerical 
scheme  developed  by  Whitfield,  Janus  and  Simpson  [WJS88],  is  a  multi-block,  finite-volume 
analysis  constructed  by  Huff,  Swafford  and  Reddy  [HSR91]  for  the  time-accurate  resolution 
of  nonlinear,  two-dimensional,  unsteady  flows  through  vibrating  cascades. 

In  this  report  numerical  results  are  presented  for  two-dimensional  unsteady  flows  through 
a  representative  high-speed  compressor  cascade,  i.e.,  the  so-called  Tenth  Standard  Cascade 
Configuration  [FS83,  FV93].  In  particular,  we  consider  unsteady  subsonic  and  transonic  flows 
driven  by  prescribed  single-degree-of-freedom,  torsional  and  bending  vibrations  of  the  blades. 
The  Euler/Navier-Stokes  results  are  compared  with  those  obtained  from  the  linearized  in¬ 
viscid  flow  analysis,  LINFLO  [Ver92]  to  help  validate  and  demonstrate  Euler/Navier-Stokes 
capabilities  for  unsteady  cascade  flows. 
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2.  Physical  Problem  and  Mathematical  Description 


2.1  Unsteady  Flow  through  a  Two-Dimensional  Cascade 

We  consider  time-dependent  flow,  at  high  Reynolds  number  (Re)  and  with  negligible 
body  forces,  of  a  perfect  gas  with  constant  specific  heats  through  a  two-dimensional  cascade, 
such  as  the  one  shown  in  Figure  1.  In  general,  the  unsteady  fluctuations  in  the  flow  can  arise 
from  one  or  more  of  the  following  sources:  blade  motions  and  aerodynamic  disturbances  at 
inlet  and  exit,  that  carry  energy  toward  the  blade  row.  In  the  absence  of  unsteady  excitation, 
the  flows  far  upstream  and  far  downstream  from  the  blade  row  are  assumed  to  be  at  most  a 
small  steady  perturbation  from  a  uniform  free  stream.  In  this  case  any  arbitrary  unsteady 
aerodynamic  motion  of  small  amplitude  can  be  represented  approximately  as  the  sum  of 
independent  entropic,  vortical,  and  acoustic  disturbances,  as  indicated  in  Figure  2.  Inlet 
entropic  and  vortical  excitations  model  the  wakes  from  upstream  blade  rows;  inlet  and  exit 
acoustic  excitations,  the  potential  flow  variations  associated  with  adjacent  blade  rows. 

In  the  present  discussion  all  physical  variables  are  dimensionless.  Lengths  have  been 
scaled  with  respect  to  blade  chord,  time  with  respect  to  the  ratio  of  blade  chord  to  up¬ 
stream  freestream  flow  speed,  density  viscosity,  and  thermal  conductivity  with  respect  to 
the  upstream  freestream  values  of  these  flow  variables,  velocity  with  respect  to  the  upstream 
freestream  flow  speed,  and  stress,  and  therefore,  pressure,  with  respect  to  the  product  of  the 
upstream  freestream  density  and  the  square  of  the  upstream  freestream  speed.  The  scalings 
for  the  remaining  variables  can  be  determined  from  the  equations  given  below,  which  have 
the  same  forms  as  their  dimensional  counterparts. 

We  will  analyze  the  unsteady  flow  in  a  blade-row  fixed  coordinate  frame  in  terms  of  the 
Cartesian  spatial  coordinates,  e.g.,  (x,  y),  where  x  and  y  measure  distances  along  and  normal 
to  the  reference  blade  chord,  and  the  time  t.  To  describe  flows  in  which  the  fluid  domain 
varies  with  time,  it  is  useful  to  consider  two  sets  of  independent  variables,  say  (x,  t)  and 
(x,t).  The  position  vector  x(x,t)  =  x  -f  lZ(x,t)  describes  the  instantaneous  location  of  a 
moving  field  or  boundary  point,  say  V,  x  refers  to  the  reference  or  steady-state  position  of 
"P,  and  'R-(x,t)  is  the  displacement  of  V  from  its  reference  position.  The  displacement  field, 
n,  is  usually  prescribed  so  that  the  solution  domain  moves  with  solid  boundaries  and  is 
stationary  far  from  the  blade  row. 

The  mean  or  steady-state  positions  of  the  blade  chord  lines  coincide  with  the  line  segments 
fj  —  ^  tan  0  -f  mG,  0  <  <  cos  0,  m  =  0,  ±1,  ±2,  . . .,  where  ^  and  fj  are  Cartesian  coordi¬ 

nates  in  the  axial  flow  and  cascade  tangential  directions,  respectively,  m  is  a  blade  number 
index,  0  is  the  cascade  stagger  angle,  and  G  =  |G1,  where  G  is  the  cascade  gap  vector  which 
is  directed  along  the  rj-axis  with  magnitude  equal  to  the  blade  spacing  (see  Figure  1).  The 
blade  motions  (see  Figure  3),  are  defined  by  Xs„  =  XB^d-'^BmC^Bm?  O5  m  =  0,  ±1,  ±2,  .... 
Here,  'J^Bm  is  a  prescribed  function  of  the  position  x,  for  x  €  B^,  and  the  time  t.  It  de¬ 
scribes  the  displacement  of  a  point  on  a  moving  blade  surface,  Hm,  relative  to  its  mean  or 
steady-state  position  on  B^,  In  the  present  study  we  are  interested  primarily  in  unsteady 
flows  excited  by  prescribed  blade  motions.  However,  entropic,  S_oo(x,t),  vortical,  ^_^(x,  t), 
and  acoustic,  p/,:poo(x,  t),  excitations,  where  the  subscripts  —00  and  -|-oo  refer  to  regions  far 
upstream  (^  ^  i-)  and  far  downstream  ((f  >  ^+)  from  the  blade  row,  respectively,  can  also 
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be  prescribed  as  functions  of  x  and  t.  These  functions  must  be  solutions  of  the  fluid  dynamic 
field  equations,  which  describe  disturbances  that  carry  energy  toward  the  blade  row. 


2.2  Fluid  Dynamic  Equations 


The  equations  that  govern  the  flow  are  determined  from  the  conservation  laws  for  mass, 
momentum  and  energy,  the  constitutive  relations  for  a  Newtonian  fluid,  and  the  thermody¬ 
namic  relations  for  a  perfect  gas.  To  describe  unsteady  flows,  it  is  convenient  to  express  the 
conservation  laws  relative  to  an  arbitrary  moving  control  volume  V{t),  which  is  bounded  by 
the  moving  control  surface  <S(x,  t),  i.e.. 


and 


^  I  pdV+  [  ~p(V  -iZ)-ndS  =  0, 
at  Jv  Js 

i-  I  pYdV  -t-  /  p[V  (8)  (V  -  it)]  ■ndS=  [  f  ndS 
dt  Jv  Js  Js 

—  [  pEjdV  +  /  PEt{Y  -it)-ndS=  f  V-f-ndS  -  [  Q  •  ndS  . 
dt  Jv  Js  Js  'ts 


(2.1) 

(2.2) 

(2.3) 


Here  p,  V  and  Et  =  E-\-V'^  j2  are  the  fluid  density,  velocity  and  specific  total  internal  energy, 
respectively,  ii  is  the  velocity  of  a  point  on  the  control  surface,  T  is  the  stress  tensor,  Q  is 
the  heat  flux  vector,  n  is  a  unit  outward  normal  vector,  and  ®  denotes  the  tensor  or  dyadic 
product  of  two  vectors. 

For  a  Newtonian  fluid  having  zero  coefficient  of  bulk  viscosity,  the  stress  tensor  (dyadic) 
is  related  to  the  fluid  (thermostatic)  pressure,  P,  and  velocity,  V,  by 


t  =  -PI  +  il  =  PI  -  pRe-^  [(2V  ■  V/3)I  -  V  0  V  -  (V  ®  V),]  ,  (2.4) 


where  I  is  the  unit  tensor;  JT  is  the  viscous  stress  tensor;  p  is  the  coefficient  of  shear  viscosity; 
Re  =  P*V*L*Ip*  is  the  flow  Reynolds  number;  the  superscript  *  refers  to  a  dimensional 
reference  value  of  a  flow  variable;  and  the  subscript  c  denotes  the  conjugate  dyadic.^  The 
heat  flux  Q  is  related  to  the  gradient  of  the  temperature  T  or  specific  enthalpy,  H  =  T,  via 
Fourier’s  law,  i.e., 

Q  =  -kPr-^Re-^YT  ,  (2.5) 

where  k  is  the  coefficient  of  thermal  conductivity,  Pr  =  p*C*/k’*  is  the  Prandtl  number  of 
the  flow,  and  C*  is  the  (dimensional)  specific  heat  of  the  fluid  at  constant  pressure.  Note 
that  because  of  the  nondimensionalization  used  in  the  present  study,  the  temperature  and 
specific  enthalpy  have  identical  numerical  values. 

The  coefficients  p  and  k  are  assumed  to  be  known  functions  of  the  temperature  T. 
For  example,  the  following  empirical  laws  relating  the  molecular  viscosity  and  the  thermal 
conductivity  to  the  temperature,  e.g.. 


r_ooj  f  +  Tc 


(2.6) 
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are  often  used  in  Navier-Stokes  calculations.  Equations  (2.6)  are  a  form  of  Sutherland’s 
Law.  Here  TLoo  =  CpTZ^I{y*^Y  is  the  nondimensional  reference  temperature,  and  Tc  is 
a  constant,  which  for  air  has  a  dimensional  value,  T^,  of  110°K  [Sch60]. 

In  addition  to  the  foregoing  equations,  two  relations  from  classical  thermodynamics  are 
required;  in  particular,  the  equation  of  state  for  a  thermally  perfect  gas  and  the  equation 
relating  the  internal  energy  and  the  temperature  or  specific  enthalpy  for  a  calorically  perfect 
gas,  i.e., 

P  =  -^ypT  and  E  =  =  (7  —  l)~^.P/p  ,  (2.7) 

respectively.  Here,  7  is  the  specific  heat  ratio  of  the  fluid  (constant  pressure  to  constant 
volume). 

The  foregoing  equations  provide  a  sufficient  set  for  determining  the  unknown  fluid- 
dynamic  variables.  For  turbulent  flows,  these  equations  describe  the  behavior  of  the  ensemble- 
or  Reynolds- averaged  values  of  the  time-dependent  flow  variables.  The  effects  of  random  tur- 
bulent  fluctuations  are  accommodated  by  incorporating  the  turbulent  correlations,  pv'  0  v' 
and  phj-y',  into  the  definition  of  the  viscous  stress  tensor  II  and  the  heat  flux  vector  Q, 
respectively.  To  obtain  closure  of  the  mean  flow  equations,  the  turbulent  correlations  are 
usually  related  to  gradients  in  the  ensemble-averaged  or  mean-flow  variables  via  algebraic 
eddy  viscosity  models. 

For  convenience  in  developing  numerical  algorithms  the  conservation  equations  are  usu¬ 
ally  written  in  column  vector  form  as 

jJ^VdV  +  J^[F,  +  G,-  -  =  0  .  (2.8) 

Here  the  state  vector  U,  and  the  flux  vectors  Fj  and  Gj,  j  =  1,2,  are  defined  by 
'  P  '  '  pVxj  '  O' 

P^Xl  pFx^Vxj  -f  Pd\j  ^XlXj 

ld=  ,  F,-=  _  ^  ^  ,  G,  =  _  .  (2.9) 

pFx2  pVx2Vxj  T  Pd2j  na:2a;j 

.  ppT  .  .  P{Et  +  P I P)Fxj  .  .  —^XjXiVxi  +  Qxj  . 
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3.  Nonlinear  Unsteady  Aerodynamic  Models 

3.1  Viscous  Flow  —  The  Reynolds-Averaged  Navier-Stokes  Equations 

Equations  that  describe  the  fluid  motion  at  the  instantaneous  position,  x  =  x  + 
of  a  moving  field  point  are  obtained  from  the  integral  conservation  laws  (2.8)  by  applying 
Reynolds’  transport  theorem  to  interchange  the  order  of  time-differentiation  and  volume 
integration,  Green’s  theorem  to  convert  surface  integrals  to  volume  integrals,  and  taking  the 
limit  of  the  resulting  volume  integrals  as  V(t)  — *•  0.  After  performing  the  necessary  algebra, 
we  find  that 

where  the  column  vectors  U,  Fj  and  Gj  are  defined  in  (2.9),  and  a  summation  over  repeated 
indices  is  implied. 

For  application  to  turbomachinery  unsteady  flows  the  field  equations  (3.1)  must  be  sup¬ 
plemented  by  near-field  conditions  at  the  vibrating  blade  surfaces  and  far-field  conditions  at 
the  cascade  inflow  and  outflow  boundaxies.  Since  transient  unsteady  aerodynamic  behavior 
is  usually  not  considered,  a  precise  knowledge  of  the  initial  state  of  the  fluid  is  not  required. 
The  no-slip  condition,  i.e.,  V  =  iZ-Bm  for  x  e  Bm  ^  where  TZ-Bm  is  prescribed,  applies  at 
blade  surfaces.  In  addition,  either  the  heat  flux  Q-n  or  the  temperature  T  must  be  prescribed 
at  such  surfaces.  The  entropy  and  vorticity  (total  temperature  and  total  pressure)  fluctua¬ 
tions  at  inlet,  and  the  incident  pressure  fluctuations  at  inlet  and  exit  must  also  be  specified. 
The  total  pressure  and  total  temperature  fluctuations  at  exit  and  the  unsteady  pressure 
response  at  inlet  and  exit  must  be  determined  as  part  of  the  time-dependent  solution. 

Computational  Coordinates 

To  enhance  the  efficiency  and  accuracy  of  a  numerical  solution  scheme  and  to  simplify 
the  implementation  of  boundary  conditions,  a  transformation  of  the  governing  equations 
from  physical  to  computational  space  is  usually  performed.  In  computational  space,  the 
solution  domain  is  a  rectangular  shape,  and  the  computational  grid  is  defined  to  be  uniform, 
orthogonal  and  stationary  for  convenience  in  defining  finite  difference  approximations. 

Thus,  we  consider  an  independent  variable  transformation  (x,  t)  (a,  r)  from  physical 
Cartesian  coordinates  x  =  (xi,a;2)  and  time  t  to  a  computational  space  described  by  the 
spatial  coordinates  cx  =  (01,0:2)  and  the  time  r.  After  performing  the  necessary  algebra,  we 
find  that  the  Navier-Stokes  equations  (3.1)  can  be  written  in  strong  conservation  form  as 

f  +  =  (3.2) 

a 

where 

0  =  =  +  ,  and  6,  =  (J'o)-' (IffG.)  . 

(3.3) 
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The  viscous  shear  stresses  [cf,  (2.4)]  and  the  heat  fluxes  [cf.  (2.5]  in  the  transformed  compu¬ 
tational  space  are  given  by 


dat  c>K,  _  Ja„dV,  • 

dxj  dak  dxi  dak  dxk  da^ 


(3.4) 


and 


-kPr-^Re-^ 


daj  dH 
dxi  daj 


(3.5) 


Values  of  the  metrics  daj/dt  and  daj/dxk  must  also  be  provided.  In  most  cases 
these  can  be  computed  numerically  using  finite-difference  approximations  for  the  derivatives 
dxk/daj,  etc.  If  the  metrics  appearing  in  Equation  (3.3)  can  be  expressed  in  terms  of  these 
derivatives,  the  numerical  computation  of  the  metrics  is  completed.  In  particular,  for  two- 
dimensional  flows  we  find  that 


dx^  da2  ’  dxr  da-,  ’ 


dx2  da2  ’  dx2  da,  ’ 


and 


da,  _ 


dx\  dx2 

dx,  dx2 

da2 

_jxa 

dx. 

dx2 

dx,  dx2 

da2  dr 

dr  da2 

’  dt 

doLi 

dr 

dr  da. 

=  1/ 


(dx,  dx2 
da,  da2 


dx,  dx2 
da2  da. 


(3.6) 


where  the  is  the  Jacobian  of  the  transformation  (x,  t)  -4  (a,r),  and  is  referred  to  as 
a  metric  Jacobian.  For  a  numerical  formulation  based  on  the  integral  form  of  the  govern¬ 
ing  equations  the  forward  metrics  are  obtained  from  the  geometry  of  computational  cell  in 
physical  space.  In  particular,  the  cell  volume  and  face  areas  are  computed  from  the  location 
of  the  grid  points  that  define  the  vertices  of  the  cell  in  physical  space. 


3.2  High  Reynolds  Number  Flow  —  The  Thin-Layer  Navier-Stokes  Equations 

The  Navier-Stokes  equations  must  be  considered  if  viscous  effects  are  expected  to  be 
important  throughout  the  fluid  domain.  However,  for  many  flows  of  practical  interest,  the 
Reynolds  number  (Re)  is  usually  sufficiently  high  so  that  the  latter  are  concentrated  in 
relatively  thin  layers  across  which  the  flow  properties  vary  rapidly  but  continuously.  Provided 
that  large  scale  flow  separations  do  not  occur,  these  layers  generally  lie  adjacent  to  the  blade 
surfaces  (boundary  layers),  downstream  of  the  blades  (wakes)  and  in  the  vicinity  of  rapid 
compressions  (shocks).  If  viscous  effects  are  negligible,  i.e.,  if  Re  — >  oo,  it  is  sufficient  to 
consider  only  the  outer  or  inviscid  flow. 

We  proceed  to  consider  high  Reynolds  number  flow  in  which  viscous  effects  are  concen¬ 
trated  in  thin  layers  lying  adjacent  to  the  blade  surfaces  and  extending  downstream  from 
the  blade  trailing  edges.  We  identify  shear  layer  surfaces  Sm  each  of  which  is  entirely  con¬ 
tained  within  a  viscous  layer,  and  postulate  that,  within  a  viscous  layer,  flow  gradients  in 
the  direction  normal  (n)  to  a  shear  layer  surface  will  be  substantially  greater  than  those  in 
the  tangential  direction.  We  can  take  advantage  of  this  feature  to  simplify  the  viscous  terms 
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that  appear  in  the  full  Reynolds-averaged  Navier-Stokes  equations.  In  particular,  we  neglect 
the  streamwise  gradients  of  the  viscous  stresses  and  the  heat  fluxes,  and  retain  the  normal 
gradients.  The  resulting  equations  are  known  as  the  thin-layer  Navier-Stokes  equations. 

The  thin-layer  equations  can  be  determined  from  the  transformed  matrix  form,  i.e.,  (3.2), 
of  the  Reynolds-averaged  Navier-Stokes  equations.  First,  we  assume  that  gradients  of  the 
viscous  flux  terms  parallel  to  a  body  surface  (i.e.,  in  the  ai -direction)  are  small  compared 
to  gradients  normal  (i.e.,  in  the  a2-direction)  to  the  surface,  and  hence  that  the  former  can 
be  neglected.  Accordingly,  we  set 

dGi  ,  . 

- —  0  (3.7) 

da-i 

In  addition,  we  assume  that  velocity  and  temperature  or  specific  enthalpy  gradients  parallel 
to  a  body  surface  are  small  compared  to  velocity  and  temperature  gradients  in  the  normal 
direction.  In  this  case  the  Cartesian  components  of  the  viscous  stress  tensor  II  and  the  heat 
flux  vector  Q  are  given  approximately  by 


.  „  _i  [ da2  dK,  da2  aK, 

^  dxj  doc2  dxi  da2  dxk  da2 


Q^.  -kPr  ^Re 


_i5q!2  dH 


dxi  da2 


respectively.  Thus,  the  viscous  flux  vector,  G2,  has  the  form  [cf.  (2.9)  and  (3.3)] 


62  =  -{jn' 


dXk^ 


dao^ 


TT 


(  Ta7a\“-1 


(3.10) 


and  it  can  be  approximated  by 


G2«g' 


\  Q0L2  dx\ 


--  (,^  da. 


(3.11) 


Ki  hh|!Zh  +  +  K,K, 

dOi2  dOi2 


where 


2 


I<i  = 


'd^ 

dxi 


+ 


'da2 

dxo 


and 


X2 


Kz  =  + 


doi2 

dxi 


dxi  da2  '  dx2  da2 

rr  da2 


(3.12) 


dx2 


The  field  equations  resulting  from  the  approximations  (3.7)  and  (3.11);  i.e.,  the  thin-layer, 
Reynolds-averaged,  Navier-Stokes  equations,  then  have  the  form 


dr 


a 


+ 


_a_ 

da' 


X  TL 

G  =  0 


(3.13) 


The  assumptions  introduced  to  arrive  at  equation  (3.13)  lead  to  a  substantial  reduction  in 
the  computer  storage  and  time  requirements  needed  to  determine  viscous  solutions. 

An  eddy  viscosity  formulation  can  be  used  in  conjunction  with  (3.11)  to  model  turbulent 
phenomena.  Thus,  an  effective  viscosity,  /igg,  and  an  effective  thermal  conductivity,  i-e.. 


)“eff  =  A  +  e  and  =  k  +  {Pr)~eH  =  +  {Prj Pr-rYe  ,  (3.14) 

are  used  in  lieu  of  fi  and  k  in  equation  (3.11).  Here  e  and  Iff  =  e/Pr^  are  the  turbulent 
eddy  viscosity  and  eddy  diffusivity,  respectively,  and  Prx  is  the  turbulent  Prandtl  number. 
The  eddy  viscosity,  e,  and  hence,  the  diffusivity  Cfi,  can  be  determined  using  an  algebraic 
turbulence  model,  e.g.,  the  two-layer,  Cebeci-Smith  [CS74]  or  Baldwin-Lomax  [BL78]  model. 


3.3  Inviscid  Flow  —  The  Euler  Equations 

In  the  inviscid  limit  {Re  oo)  the  thicknesses  of  the  viscous  layers  become  zero,  and  they 
can  be  modeled  as  surfaces  across  which  the  flow  variables  are  discontinuous.  In  particular, 
boundary  layers  and  wakes  become  vortex  sheets  that  support  a  discontinuity  in  tangential 
velocity,  and  shocks  become  surfaces  that  support  a  discontinuity  in  normal  velocity.  It  is 
usually  assumed  that  the  boundary  layers  remain  attached  to  the  blade  surfaces  and  support 
jumps  in  velocity  from  zero  at  the  “walls”  to  inviscid  values  at  their  “edges.”  Furthermore, 
the  vortex-sheet  unsteady  wakes  emanate  from  the  blade  trailing  edges  and  extend  infinitely 
far  downstream.  The  integral  forms  of  the  conservation  laws  are  required  to  describe  the  flow 
over  the  entire  domain  of  interest.  These  provide  differential  equations  in  regions  where  the 
flow  variables  are  continuously  differentiable  and  “jump”  conditions  at  the  surfaces  across 
which  the  flow  variables  are  discontinuous. 

The  field  equations  that  govern  continuous,  inviscid,  fluid  motion  (i.e.,  the  Euler  equa¬ 
tions)  are  obtained  from  equations  (2.1)-(2.3)  by  setting  T  =  —PI  and  Q  =  0,  or  from 

equations  (3.2)  and  (3.13)  by  setting  G  and  G  =0,  respectively.  In  principle,  the  in¬ 
viscid  field  equations  must  be  supplemented  by  jump  conditions  that  apply  across  moving 
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boundary  layers,  Brm  wakes,  W^,  and  shocks,  Shm-  However,  the  usual  practice  in  numerical 
computations  is  to  solve  the  field  equations  over  the  entire  fluid  domain,  thereby  “capturing” 
wake  and  shock  phenomena.  The  inviscid  flow  is  then  determined  subject  to  flow  tangency 
conditions,  i.e., 

=  0  for  xeBm.  (3.15) 

at  the  moving  blade  surfaces,  and  far-field  conditions  at  the  inflow  and  outflow  boundaries 
of  the  computational  domain.  The  far-field  conditions  used  in  the  inviscid  approximation 
are  the  same  as  those  indicated  previously  for  Navier-Stokes  simulations. 
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4.  The  NPHASE  Analysis 


A  numerical  model,  based  on  Roe’s  [RoeSl]  approximate  Riemann  solver,  for  solving  non¬ 
linear  unsteady  flow  problems  was  developed  by  Whitfield,  Janus  and  Simpson  [WJS88].  This 
model  was  subsequently  extended  for  turbomachinery  applications  and  implemented  into  an 
unsteady  flow  code,  called  NPHASE,  by  Huff,  Swafford  and  Reddy  [HSR91].  NPHASE  is 
an  implicit,  multi-block,  cell- centered,  finite-volume  code,  that  can  be  used  for  predicting 
two-dimensional,  nonlinear,  viscid  (at  high  Reynolds  number)  and  inviscid  unsteady  flows 
through  vibrating  cascades.  A  detailed  description  of  the  current  version,  i.e.  version  2.0,  of 
the  NPHASE  analysis  can  be  found  in  [SLH+94]  and  the  references  cited  therein. 

The  computational  mesh  used  in  NPHASE  is  a  sheared  H-mesh,  typically  generated  using 
the  IGB  grid  generation  package  of  [BH92].  This  structured  mesh  defines  a  curvilinear  coor¬ 
dinate  system,  in  which  coordinate  curves  lie  along  the  boundaries  of  the  physical  domain. 
A  time- dependent  coordinate  transformation,  (x,  t)  from  the  physical  domain,  in 

which  the  grid  deforms  with  the  blade  motion,  to  a  computational  domain,  in  which  the 
grid  is  stationary,  uniform,  and  orthogonal,  is  applied  to  simplify  the  implementation  of 
numerical  differencing  and  flow  boundary  conditions.  Recall  that  x  =  x  -1-  7?.(x,  t),  where  TZ 
defines  the  displacement  of  the  moving  field  points  from  their  reference  or  mean  positions. 

The  time-dependent  displacement  field,  7Z{5c,  t),  is  prescribed  in  a  manner  that  facilitates 
implementation  of  the  flow  boundary  conditions.  It  is  determined  here  as  a  solution  of 
Laplace’s  equation,  =  0;  subject  to  Dirichlet  boundary  conditions  at  the  blade  surfaces, 

i.e.,  71  =  TZ-Bm  for  X  G  Bm,  and  in  the  far  field,  7?.  =  0  for  ^  and,  for  the  blade  motions 
considered  herein  [cf.  (5.1)-(5.3)],  phase-lagged,  blade-to-blade,  periodicity  conditions;  i.e., 
-h  mG,t  —  mcrlu})  =  7?.(x,  t),  upstream  and  downstream  of  the  blade  row. 

4.1  Finite  Volume  Approximation 

For  a  finite  volume  discretization,  based  on^  the  integral  forms  of  the  governing  field 
equations,  the  time-dependent  geometric  properties  of  the  mesh  cells  in  physical  space  are 
required.  These  include  the  cell  volume,  t?,  the  volume  swept  out  per  unit  time  by  the 
constant  aj  face  as  the  cell  interface  moves,  ■i9j,  and  the  area  of  the  constant  qj  face  projected 
in  the  Xk  direction,  Ajk-  These  geometric  properties  of  a  cell  are  determined  from  the 
instantaneous  locations  of  the  cell  vertices  in  physical  space. 

We  can  write  the  finite  volume  spatial  discretization  of  equation  (3.13)  as 

=  -  82G^  =  -i.  (4.1) 

a 


dr 


where 

U  =  ,  Fj  =  -f-  AjfcFfc  ,  and  G  =  {r‘^)-^A2kGk  •  (4.2) 

Here,  the  vector  U  represents  an  average  over  the  cell  volume;  the  vector  Fj  is  the  inviscid  flux 

»  TZ,  . 

across  a  constant  aj  cell  face;  G  is  the  viscous  flux  across  a  constant  a2  cell  face;  the  symbol 
6j  denotes  the  difference  across  adjacent  cell  interfaces,  i.e.  6j{  )  =  (  )j+i/2  —  (  )j-i/2'i 
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and  R  is  referred  to  as  the  residual.  The  repeated  j  index  in  (4.1)  implies  summation 
over  all  computational  coordinate  directions,  so  that  the  term  SjFj  is  the  net  inviscid  flux 
through  the  cell.  The  grid  deforms  in  physical  space  as  the  blades  move.  Therefore,  the  cell 
geometry  terms,  dj  and  Ajk,  are  time  dependent.  A  comparison  of  equations  (3.3)  and 
(4.2)  indicates  that  these  geometric  terms  can  be  interpreted  in  terms  of  the  metrics  of  the 
(x,t)  — >  (q:,t)  coordinate  transformation.  In  particular,  t?  =  (J^“)~^,  d  =  {J^°')~^daj/dt 
and  Ajk  =  {J^^^Y^dajj dxk- 

In  the  NPHASE  analysis,  the  fluxes  at  the  cell  interfaces  are  computed  from  the  values 
of  the  state  vector  on  either  side  of  the  interface  using  the  approximate  one-dimensional 
Riemann  solver  developed  by  Roe  [RoeSl].  The  interfacial  flux  is  computed  by  multiplying 
the  difference  in  the  state  vector  across  the  interface  by  a  flux  Jacobian  matrix  representing 
the  local  interface  conditions.  The  eigenvalues  of  this  flux  Jacobian  matrix  are  used  to 
determine  which  characteristic  modes  axe  included,  thus  controlling  the  direction  of  spatial 
differencing.  This  technique  is  known  as  flux  difference  splitting  and  results  in  a  first-order 
accurate  spatial  differencing  scheme.  Higher  order  accuracy  is  obtained  by  using  a  corrective 
flux,  which  is  limited  by  a  TVD  scheme  to  control  dispersive  errors. 

The  NPHASE  time  differencing  is  based  on  a  second-order,  implicit,  three-point,  back¬ 
ward,  difference  approximation  of  the  time  derivative.  Thus,  if  the  superscript  n  refers  to 
the  current  or  nth  time  level,  we  can  write 


■z  n  z  n— 1 

3AU  -AU 
2At 


-R 


n+1 


(4.3) 


Z  Tl  z  Tl-rx  ^  Tl 

where  AU  =  U  —  U  ,  and  the  differences  on  the  left-hand  side  involve  both  the  state 
vector  and  the  cell  volume.  Separating  the  time  dependence  of  the  state  vector  and  the  cell 
volume  results  in 

—  «_}_1  ^  ^ 

0AU"  +  R  =  ??"-^AU""V2Ar  -  U”(3i?”+^  -  4??"  -f  t?"~^)/2Ar  =  ,  (4.4) 

where  ©  =  3i?”"'’^/2At. 

The  nonlinear  equation  (4.4)  is  solved  at  each  time  step  using  a  Newton  iteration  proce¬ 
dure  in  which  the  iteration  formula  is  given  by 


01  + 


5R 


dV 


AW  =  0(UP-\-U’^) 


-  p— 1  ,1 

R  + 


(4.5) 


Here  p  —  1,2, .. .  is  the  Newton  iteration  index,  U°  =  U”,  AU^  =  W  —  W  is  the 

Newton  update  to  the  state  vector,  and 


dR 

dlJ 


AW  =  Sj 

Up-1 


diJ 


tjp-i 


(4.6) 


Once  the  Newton  iteration  converges  AU^  =  0,  and  (4.4)  is  satisfied.  Note 

that  in  the  current  version  of  NPHASE  [SLH+94]  the  viscous  flux  term  is  included  in  the 

~  p— 1 

right-hand  side  residual  R  of  equation  (4.5),  but  ignored  in  constructing  the  left-hand 
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side,  implicit,  operator  5R/5tJ  yp_j-  This  approach  introduces  errors  into  the  left-hand  side 
operator,  but  these  become  negligible  as  the  Newton  iterations  converge. 

The  inviscid  flux  terms  on  the  right-hand  side  of  (4.5)  are  evaluated  using  flux  difference 
splitting  corrected  for  higher  order  spatial  accuracy.  Those  on  the  left-hand  side  [cf.  (4.6)] 
are  evaluated  using  flux  vector  splitting  for  convenience  in  defining  the  approximate  factor¬ 
ization,  used  to  facilitate  the  solution.  The  coupling  between  adjacent  cells  introduced  by 
the  flux  terms  on  the  left-hand  side  of  (4.5)  makes  it  expensive  to  solve  this  equation.  So  the 
flux  evaluation  is  approximately  factored  into  the  product  of  positive  and  negative  operators. 
To  reduce  the  error  introduced  by  this  approximate  factorization,  symmetric,  Gauss-Seidel, 
subiterations  [WT91]  are  used  as  part  of  the  time  stepping  procedure.  The  subiteration 
procedure  employs  an  LU  decomposition  of  the  Newton  iteration  matrix,  with  forward  and 
backward  substitution. 

The  field  equation  (4.1)  must  be  solved  subject  to  the  appropriate  surface  conditions  at 
the  blade  boundaries,  and  far-field  conditions  at  the  inflow  and  outflow  boundaries  of  the 
computational  domain.  The  far-field  conditions  should  permit  the  prescription  of  external 
unsteady  aerodynamic  excitations,  and  allow  unsteady  disturbance  waves  coming  from  within 
the  solution  domain  to  pass  through  the  inflow  and  outflow  boundaries  without  reflection. 
The  flow  tangency  condition  used  in  the  inviscid  version  of  the  NPHASE  analysis  is  based  on 
a  two  phantom-cell,  pressure-symmetry,  implementation  to  lower  the  generation  of  spurious 
numerical  entropy  and  vorticity  at  a  blade  surface.  In  an  earlier  version  [HSR91]  of  the 
NPHASE  analysis,  the  computational  mesh  was  stretched  in  the  axial  direction  to  dissipate 
outgoing  waves.  Recently,  however,  far-field  conditions  based  on  one-dimensional  charac¬ 
teristic  theory  and  on  approximate  two-dimensional  characteristic  theory  [Gil90],  have  been 
implemented.  Although,  these  conditions  allow  for  the  prescription  of  incoming  unsteady 
excitations,  they  result  in  reflections  of  outgoing  disturbances,  if  the  latter  have  significant 
variation  in  the  circumferential  direction.  As  a  result,  the  stretched  mesh  capability  was 
applied  in  the  NPHASE,  blade  vibration,  calculations  reported  herein. 


4.2  Flux  Evaluations  and  Approximate  Factorization 


The  flux  splitting  used  in  the  NPHASE  analysis  is  based  on  a  similarity  transformation 
and  an  eigenvalue  decomposition,  i.e.. 


df  _  / dfy  / \ 


t(A^  -h  A  )T-'  , 


(4.7) 


of  the  flux  Jacobian  matrix,  dF/dU,  into  matrices  that  account  for  right  (-f)  and  left  (— ) 
traveling  disturbances.  The  matrices  T  and  contain  the  right  and  left  eigenvectors, 
respectively,  of  dF/dlJ,  and  A  and  A  are  diagonal  matrices  containing  the  positive  and 
negative  eigenvalues. 

Flux  vector  splitting  is  used  to  approximately  factor  the  left-hand  side  of  (4.5).  This 
approximation  has  the  form  [SW81] 


dF  ^  ~  dF 

—  ^  AXJ -|-  — ~ 

dV  dV 


(4.8) 
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where  the  subscripts,  /,  L,  and  R  refer  to  a  cell  interface  and  the  cell  volumes  to  the  “left” 
and  “right”  of  this  interface,  respectively.  It  results  in  first  order  spatial  accuracy,  but  is 
only  used  in  the  approximate  factorization,  and  therefore,  does  not  appear  in  the  converged 
final  solution. 

It  follows  from  (4.8)  that 


(4.9) 


The  approximate  factorization  equation  for  the  Newton  iteration  procedure  is  therefore 

5aup  -  Au^_i  +  m“+i A%i  =  -0(up-^  -  u")  -  ,  (4.10) 

where  j  is  a  grid  point  index  corresponding  to  the  jth.  computational  coordinate  direction 
and  p  is  the  Newton  iteration  index.  The  D  and  M  matrices  are  evaluated  based  on  the  state 
vector  at  the  previous  Newton  iteration.  The  D  matrix  contains  the  diagonal  elements 
of  the  iteration  matrix,  and  the  M'*'  and  M~  matrices  contain  the  off-diagonal  elements  in 
the  negative  and  positive  computational  coordinate  directions,  respectively.  These  matrices 
are  given  by 

: ). 

Up-1/ 

(4.11) 

Up-1 

Equation  (4.10)  is  solved  for  AU^  using  a  symmetric  Gauss-Seidel  subiteration  procedure. 
After  introducing  I  as  the  Gauss-Seidel  iteration  index,  we  find  that  the  subiteration  formula 
is  given  by 

DAU^’*  -  M^t.1  AU^fi  + 

DAU^-'  +  m”+i AU^ii  -  I^^tiAU^l*! 

The  first  subiteration  is  over  positive  grid  indices;  the  second,  over  negative  grid  indices.  The 
subiteration  procedure  is  thus  an  LU  decomposition  of  the  Newton  iteration  matrix,  with 
forward  and  backward  substitution.  Once  the  Gauss-Seidel  subiteration  procedure  converges, 
equation  (4.5)  is  satisfied. 

As  the  solutions  converge,  any  errors  introduced  by  the  Newton  iteration  or  the  ap¬ 
proximate  factorization  vanish.  Only  the  errors  in  the  residual  calculation  of  equation  (4.1) 
remain.  The  inviscid  flux  terms  in  the  residual  are  calculated  using  flux  difference  splitting 


=  -0(UP-i  -  U”)  -  R^"  -h 

(4.12) 

=  -©(tJP-^  -  U”)  -  r'’  V 


=  0I+E 

i=i 


dFi 


(§  AW 

\  Up-1 


4  aw)  4^' 


JJp-l 


with  a  corrective  flux,  as  given  by  (4.16),  added  to  obtain  higher  order  spatial  accuracy.  The 
flux-difference  splitting  approximation  for  the  inviscid  flux  terms  has  the  form  [RoeSl] 


F,  =  F(Ui)  +  ^  (Uh  -  Ut)  .  (4.13) 

Ur.. 

where  the  intermediate  state  vector,  Uroc)  is  defined  using  the  relations 

-  j  77  \lfL^T,L  +  \I}rEt,r 

PRoe  —  V  PlPR  •>  *  Roe  —  ?  and  £'J’,Roe  —  r£-  — 

ypL'^  \J  Pr  V  Pl  +  V  Pr 

(4.14) 

Once  the  inviscid  fluxes  have  been  computed,  they  are  spatially  differenced  to  compute 
the  net  flux  through  the  control  volume.  The  difference  expression  is 


F 


(4.15) 


where  j  corresponds  to  the  computational  coordinate  direction,  and  the  fractional  grid  indices 
indicate  cell  interfaces.  The  difference  approximation  (4.15)  is  first-order  accurate  if  the 
flux  is  based  only  on  information  from  adjacent  cells,  but  higher  order  spatial  accuracy 
can  be  achieved  by  adding  a  corrective  flux.  The  corrective  flux  brings  in  information 
from  additional  neighboring  cells,  but  requires  the  use  of  a  flux  limiter  to  control  dispersive 
errors.  Various  flux  limiting  schemes  exist,  and  the  choice  between  them  is  not  clear.  The 
NPHASE  code  currently  supports  min-mod  [OC84],  Superbee  [Roe85],  and  Van  Leer  flux 
limiters  [VL74],  with  Van  Leer  being  preferred.  The  limiters  are  activated  by  changes  of  sign 
in  the  jumps  in  characteristic  variables  at  adjacent  interfaces,  such  as  occur  at  shocks  and 
at  stagnation  points. 

With  the  addition  of  a  corrective  flux,  the  flux  vector  has  the  form 


Fi  =  F(Ur)  + 


(Uj,  -  Ui.)  +  jT(UR„,)l£r+  -  » 


(4.16) 


where  the  limited  jumps  in  the  characteristic  variables  are  given  by 

cr+  =  t(&f_^,  erf)  and  i-'  =  L(i-J  . 


(4.17) 


Here  /  —  1  and  I  +  1  refer  to  the  cell  interfaces  to  the  “left”  and  “right”  respectively,  of 
the  interface,  i.e.,  I,  under  consideration.  The  jumps  in  the  characteristic  variables  at  the 
I  -t-  M,  M  =  —1, 0, 1,  cell  interface  are  given  by 

=  A*(URoe)T-i(URoe)  (Ur+m  -  Vl+m)  (4.18) 


^I+M  —  ^  (URoe)T  ^(IJRoe)  (^Ur+m  —  Ur+a/J  (4-18) 

where  R  +  M  and  L  +  M  are  the  cell  centers  to  the  right  and  left,  respectively  of  the  interface 
I  -|-  M.  The  matrices  T  and  A  in  (4.18)  are  computed  based  on  Roe  averaged  variables,  and 
the  function  L  in  (4.17)  is  a  generic  limiting  function.  For  example,  the  Van  Leer  limiter  is 
defined  by 


^n  +  ^n 


(4.19) 


where  the  subscript  n  refers  to  the  nth  component  of  the  column  vector. 
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4.3  Viscous  Terms 


The  viscous  flux  terms  appearing  in  equation  (4.5)  are  evaluated  at  the  cell  interfaces,  I, 
at  each  step  of  the  Newton  iteration  procedure.  The  viscous  fluxes  are  based  on  the  values  of 
state  vector,  which  are  known  from  the  previous  Newton  iteration,  in  the  cell  volumes 

to  the  left,  L,  and  right,  R,  of  the  interface,  i.e., 

_  —  'j'j^ 

G  =G  (Ur'.UD-  (4-20) 

I 

The  individual  terms  in  the  viscous  flux  vector  [cf.  (3.11)]  are  considered  separately.  The  grid 
metrics,  dcx^ldxj^  are  known  at  the  current  time  level,  and  the  velocity  derivatives  dVxJdct^ 
are  evaluated  at  the  cell  interface  using  a  central  difference  approximation,  i.e., 


— —  =  Vt 

da2  r  ^  R 


(4,21) 


The  velocity  components,  and  the  metric  Jacobian,  must  also  be  evaluated  at  cell 
interfaces.  These  quantities  are  determined  by  averaging  their  values  at  the  adjacent  cell 
volumes  and  are  given  by 

Once  the  viscous  fluxes  have  been  computed,  they  are  spatially  differenced  according  to 


~TL  -TL  ~TL 

62G  =  -  Gj_i 


(4.23) 


to  compute  the  net  flux  through  the  control  volume. 


Turbulence  Model 

An  eddy  viscosity  formulation  has  been  incorporated  into  the  NPHASE  analysis  to  model 
turbulent  phenomena.  The  eddy  viscosity,  e,  and  hence,  the  diffusity,  is  deter¬ 

mined  using  the  Baldwin  Lomax  model  [BL78]  for  blade  boundary  layers  and  the  Thomas 
model  [Tho79]  for  wakes.  Thus,  at  blade  surfaces  we  set 

f.  Y  <Y 

(4.24) 

^outer  ^  ^  pCTOSsover 

where  the  subscripts  refer  to  the  inner  and  outer  turbulent  regions  and  to  the  interface 
between  these  regions,  and  Y  and  i^rossover  are  normal  distances  from  a  blade  surface.  Y  is 
the  normal  distance  to  a  cell  center  and  L^rossover  is  the  smallest  normal  distance  at  which 
the  inner  and  outer  values  of  the  turbulent  eddy  viscosity  are  equal. 

The  turbulent  eddy  viscosity  in  the  inner  region  of  a  boundary  layer  is  determined  using 
the  Prandtl-Van  Driest  formulation  and  is  given  by 

(4.25) 
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where 

i  =  XY  \l-exp  (-r+/A+)] 

is  a  characteristic  length  scale  and 

dx2  dxi 


(4.26) 

(4.27) 


is  the  magnitude  of  the  vorticity.  In  (4.26),  X  =  0.4  is  von  Karman’s  constant,  =  26,  and 
iP'wi  where  the  subscript,  to,  indicates  that  the  quantities  are  evaluated  at 
the  wall  or  blade  surface.  The  wall  shear  stress  ,  is  given  by 


(4.28) 


where  I4i  is  the  velocity  component  tangential  to  a  blade  surface,  i.e.,  in  the  cci-coordinate 
direction,  and  the  partial  derivative  is  taken  normal  to  the  blade  surface. 

The  eddy  viscosity  in  the  outer  region  of  a  boundary  layer  is  given  by 

^outer  ~  ^cppF wake^^Kleb  >  (4.29) 

where  the  Clauser  constant,  /f',  and  the  constant  Cep  are  set  equal  to  0.0168  and  1.6,  respec¬ 
tively, 

=  min  [y„,axFmax,  %^(Knax  -  Knin)"]  ,  (4.30) 

^  max  -1 

the  Klebanoff  intermittency  factor,  FKieb,  is  given  by 

i=Kieb=  l  +  5.5(^^)  ,  (4.31) 

and  the  constant  CKieb  is  set  equal  to  0.3.  In  equation  (4.30)  Fmax  and  are  determined 
by  the  function 

F{Y)  -  y|<S|[l  -  exp(-y+M+)]  .  (4.32) 

where  jPmax  =  T’(i7iiax)  is  the  maximum  value  of  F  at  the  streamwise  station  {cxi  =  constant) 
under  consideration,  and  the  constant  C^k  is  set  equal  to  0.25. 

The  wake  eddy  viscosity  is  governed  by  the  Thomas  model  [Tho79],  in  which  the  eddy 
viscosity  is  defined  as  in  (4.25),  but  the  characteristic  length  scale,  £,  is  given  by 

i  =  4(ymax  -  Knin)/|5Uax  •  (4.33) 

The  constant  io  in  (4.33)  is  set  equal  to  0.09,  t4iax  and  t4jn  are  the  maximum  and  minimum 
values,  respectively,  of  the  flow  speed,  and  l^lmax  is  the  maximum  value  of  the  vorticity,  at 
the  streamwise  (oi  =  constant)  station  under  consideration. 
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5.  Unsteady  Excitation  and  Aerodynamic  Response 


5.1  Blade  Motion 


We  will  restrict  our  consideration  to  unsteady  flows  in  which  each  two-dimensional  blade 
undergoes  a  rigid-body  motion  or  the  form 

- 1]  •  Rp  =  hm(0  +  (ez  X  Rp)smam{t)  -  Rp[l  -  cosQ;,„(i)]  ,  (5.1) 

where 

cos  dm (t)  -  sin  am{t) 

0^{t)=  ,  (5.2) 

_  sin 6:^(0  cosam(t) 

hm(t)  and  &Tn{t)  =  &Tn{t)ez  are  the  translational  and  rotational  displacement  vectors,  re¬ 
spectively,  for  the  mth  blade,  Rp  is  the  distance  from  the  steady-state  position  of  the  mth 
blade  pitching  axis  at  Xp^  =  xp  -|-  mG  to  the  point  XB,n  =  xb  +  mG,  and  is  a  unit  vector 
that  points  out  from  the  page. 

In  addition,  we  assume  that  the  translational  and  rotational  displacement  vectors  are  of 
the  form 

hm  =  Re{hexj>[i{ut  +  mcr)]}  and  =  Re  {a  exp  [si  (wt  -|-  mcr)]}  ,  (5.3) 


where  h  and  o:  are  the  complex  amplitudes  of  the  reference  (m  =  0)  blade  translational  and 
rotational  displacements,  respectively,  u  is  the  frequency  of  the  blade  motion,  a  is  the  phase 
angle  between  the  motions  of  adjacent  blades,  and  Re{  }  denotes  the  real  part  of  {  }. 

In  general,  the  components  of  the  displacement-amplitude  vectors  h  and  a  are  complex  to 
permit  phase  differences  between  the  translational  and  the  rotational  modes  of  blade  motion. 

After  expanding  the  trigonometric  functions  in  (5.1)  into  power  series  in  dm,  we  find  that 
the  rigid  body  motion  of  the  mth  blade  can  be  expressed  in  the  form 


(-irag* 

(2n  +  l)! 


+  RpE 


Here 

^Bm  =  hm  +  am  X  Rp  =  Re{(h  -f  d  x  Rp)  exp[i(u;t  +  ma)]}  (5.5) 

is  the  approximate  representation  of  a  rigid-body  blade  motion  that  is  used  in  linearized 
unsteady  aerodynamic  analyses. 


5.2  Unsteady  Aerodynamic  Response  at  a  Moving  Blade  Surface 

A  fluid-dynamic  boundary  value  problem,  e.g.,  one  of  those  posed  in  §  3,  must  be  solved 
numerically  to  determine  the  unsteady  aerodynamic  response  information  needed  for  blade 
row  aeroelastic  and  aeroacoustic  applications.  In  particular,  for  aeroelastic  applications, 
information  at  the  moving  blade  surfaces,  e.g.,  surface  pressures  and  global  unsteady  airloads, 
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is  required.  The  pressure  acting  at  the  moving  reference  blade  surface  B  can  be  expressed 
in  terms  of  the  independent  variables  xg  and  t,  or  Xb  and  t,  i.e., 

Pb  =  =  Pb[xb  +  72.b(xb,05 =  -Pe(xB,t)  .  (5.6) 

It  is  usually  assumed  that  a  periodic  unsteady  excitation  with  temporal  period  27r/i(;  will 
produce  an  unsteady  pressure  response  at  the  moving  blade  surface  B  that  is  also  periodic 
in  time,  with  period  2Tr/(jj.  In  this  case  we  can  write 

oo  oo 

P^(xB,i)=  P„(xB)exp(ma;t)  =  Po(xb)  +  X]  Pe{p„(xB)exp(za;i)}  ,  (5.7) 

n=-oo  n=l 

where  the  Fourier  coefficients,  Pn,  are  given  by 

Pn(xB)  =  (Stt)  ^  PB{xB,t)exTp{-inLot)d{iot),  n  =  0,  ±1,  ±2, . . .  .  (5.8) 

Here,  Pq  is  the  temporal  mean  value  of  the  unsteady  pressure  and  =  2P„,  n  =  1,2,...  is 
the  complex  amplitude  of  the  nth-harmonic  (n  >  1)  component  of  the  unsteady  pressure. 

The  global  unsteady  air  loads  are  determined  by  performing  the  appropriate  integra¬ 
tions,  involving  the  unsteady  pressure,  over  a  moving  blade  surface.  For  example  the  time- 
dependent  force  and  moment  acting  on  the  reference  blade  are  given  by 


F  =  —  ^  pBllBdT B 


(5.9) 


Pb^b  X  R7> 


nsdrs 


(5.10) 


where  ng  is  a  unit  normal  vector  pointing  out  from  the  moving  blade  surface,  B,  and  is 
a  differential  distance  along  this  surface,  taken  as  positive  in  the  counterclockwise  direction. 
In  equation  (5.10)  the  moment  is  taken  about  the  moving  pitching  axis  located  at  x  =  x-p  = 
xp  -f  h  and  Rp  =  xb  —  xp. 


Aerodynamic  Work  per  Cycle 


Two  other  surface  response  parameters,  that  are  useful  in  blade  aeroelastic  investiga¬ 
tions,  are  the  aerodynamic  work  per  cycle,  Wc,  and  the  local  work  per  cycle  or  pressure- 
displacement  function,  wc{x-b)  [Ver89a].  The  aerodynamic  work  per  cycle,  Wc,  represents 
the  work  done  by  the  fluid  on  a  given  blade  over  one  cycle  of  its  motion.  In  particular,  a 
prescribed,  blade  motion  can  be  classified  as  stable,  neutrally  stable  or  unstable  depending 
upon  whether  the  aerodynamic  work  per  cycle  is  less  than,  equal  to,  or  greater  than  zero, 
respectively.  The  pressure-displacement  function  describes  the  distribution  of  the  work  per 
cycle  over  a  blade  surface.  By  definition 

I  dw  1  pij>+2ir  j•(f>+2^^  />  .  OTZ-u 

=  =  (5.11) 

where  W  is  the  instantaneous  work  done  by  the  airstream  on  the  reference  blade,  and  dW  jdt 
is  the  rate  at  which  this  work  is  done. 
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If  the  reference  blade  displacement,  TZ-b,  and  surface  pressure,  Pg,  are  expressed  as 
functions  of  the  independent  variables  xb  and  t,  as  in  (5.4)  and  (5.6),  then  an  integral  over 
the  moving  blade  surface  B  can  be  referred  to  the  mean  or  steady-state  position,  B,  of  this 
surface.  In  this  case,  the  order  of  the  integrations  used  to  determine  the  work  per  cycle  can 
be  interchanged.  For  a  rigid  blade  motion  dr^  =  dfs,  and 

Wc  =  ^  ^  wc{r^)dfB  .  (5.12) 


The  function 


r'^+2^  5  u 

Pb—^  ■  nBd{ujt} 


is  the  local  work  per  cycle  or  pressure-displacement  function.  The  unit  normals  and  hb 
at  the  moving  {B)  and  steady-state  (B)  positions,  respectively,  of  the  reference  blade  surface, 
are  related  by 


-tb  + 


_  ,  OIZb 

hb  +  -XI—  X  e, 

otb 


(6.14) 


Here  ng  and  Hb  point  outwaxd  from  the  blade  surface  and  the  unit  tangent  vector,  tb,  points 
in  the  counterclockwise  direction. 


Small- Amplitude  Blade  Motions 

If  we  restrict  our  consideration  to  rigid-body  blade  displacements  of  small-amplitude,  we 
can  set 

”^B  =  Jb  ~  q:^R-p/2  ,  (5.15) 

ng  =  hB(l  -  tt^/2)  +  octb  +  . . .  ,  (5.16) 

and 

•  hb  =  dfB  ‘  Tb  +  (fe^  ddR?)  •  db  +  •  •  •  ,  (5-17) 

where  the  dots  refer  to  terms  that  are  of  0(|72.bP)  ot  higher.  If  h  and  dc  are  harmonic  in 
time  with  frequency  a;  as  in  equation  (5.3),  then 

•  ns  =  ^  y)  [Pm  exp{imut)  -  R*^  exp(-zma;t)]  -h  .  •  •  ,  (5.18) 

^  m=0 

where  Rq  —  Q;*rB  •  "^6/2,  R\  ■=  Vb  •  hb,  and  R2  =  (Q:rB  •  tb  —  •  hB)/2.  After  substi¬ 

tuting  (5.18)  and  the  Fourier  series  expansion  for  the  surface  pressure,  equation  (5.7),  into 
equation  (5.13)  and  carrying  out  the  integrations  with  respect  to  ut,  we  find  that 

Wc  =  —TrIm[Poar^  •  tb  +  PiTg  •  hb  +  p2(Q:*rB  •  -  Q;*^Rp  •  hB)/2]  +  . . .  .  (5.19) 

If,  as  in  linearized  unsteady  aerodynamic  calculations,  we  assume  that  the  steady  pressure 
is  of  0(1)  and  the  nth  harmonic  of  the  unsteady  pressure  is  of  0(|7?.b|”),  then 


Wc  =  — 7r/m{PoQ:rB  •  tb  +  Piv’^  ■  hb}  +  0{\'Rb\‘^)  ■ 


(5.20) 


Thus,  the  higher  harmonic  content  (n  >  2)  of  the  unsteady  flow  would  have  a  negligible 
impact  on  the  energy  transfer  from  the  airstream  to  the  rigid-body  blade  motion.  Even  if 
the  unsteady  components,  n  =  1,2, 3 . . .  of  the  fluid  pressure  are  all  of  the  same  order, 
i.e.,  C>(|7?.b|),  then  to  leading  order,  i.e.,  the  pressure-displacement  function  would 

depend  only  upon  the  steady  and  the  first  harmonic  components  of  the  fluid  pressure.  Thus, 
for  harmonic,  rigid  body,  blade  motions  of  sufficiently  small  amplitude,  the  higher  harmonic 
content  of  the  unsteady  pressure  response  should  have  a  negligible  impact  on  the  transfer  of 
energy  from  the  fluid  to  the  structural  motion.  This  analytical  result  can  be  used  to  guide 
the  interpretation  and  validation  of  the  unsteady  aerodynamic  response  predictions  that  are 
determined  using  the  time-accurate,  nonlinear,  analysis  NPHASE. 
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6.  Numerical  Results 


Unsteady  aerodynanaic  response  predictions  will  be  presented  for  two-dimensional,  invis- 
cid  and  viscid  (at  Re  =  10®),  flows  through  a  compressor-type  (M+oo  <  M-oo)  cascade  to 
demonstrate  and  validate  the  current  capabilities  of  the  NPHASE  analysis  and  code.  The 
cascade,  known  as  the  Tenth  (10th)  Standard  Configuration  [FS83,  FV93],  consists  of  mod¬ 
ified  NACA  5506  airfoils,  i.e.,  the  airfoil  thickness  distribution  has  been  modified  so  that 
the  blades  close  in  wedge-shaped  trailing  edges.  It  has  a  stagger  angle,  0,  of  45  deg  and  a 
gap/chord  ratio,  G,  of  unity,  and  operates  at  uniform  subsonic  inflow  conditions.  We  will 
consider  unsteady  flows  through  the  10th  Standard  Cascade  operating  at  an  inlet  freestream 
Mach  number  M_oo  and  flow  angle,  O_oo,  of  0.7  and  55  deg  respectively,  and  at  an  inlet 
freestream  Mach  number  and  flow  angle  of  0.8  and  58  deg.  In  the  first  case  the  steady  back¬ 
ground  flow  is  entirely  subsonic;  in  the  second,  it  is  transonic  with  a  normal  shock  emanating 
from  the  suction  surface  of  each  blade. 

In  particular,  we  will  consider  unsteady  flows  excited  by  prescribed  single-degree-of- 
freedom  (SDOF),  harmonic,  blade  motions.  The  latter  occur  at  unit  frequency  and  with 
constant  phase  angle  a  between  the  motions  of  adjacent  blades.  The  motions  considered  are 
pure  translations  normal  to  the  blade  chords  (bending)  and  pure  rotations  (torsion)  about 
axes  at  the  blade  midchords.  The  blade  motions  are  termed  subresonant  if  all  acoustic 
response  disturbances  attenuate  with  increasing  axial  distance  from  the  blade  row,  superres¬ 
onant  (m,  n)  if  m  and  n  such  disturbances  persist  in  the  far  upstream  and  far  downstream 
regions,  respectively,  and  carry  energy  away  from  the  blade  row,  and  resonant  if  at  least  one 
acoustic  response  disturbance  persists  in  either  the  far  upstream  or  far  downstream  regions 
of  the  flow  and  carries  energy  along  the  blade  row  [Ver89b]. 

In  addition  to  the  NPHASE  results,  for  purposes  of  comparison,  we  will  also  present 
linear  response  predictions  as  determined  using  the  LINFLO  analysis.  In  the  latter,  the 
unsteady  flow  is  regarded  as  a  small  inviscid  perturbation  of  a  potential  steady  background 
flow.  The  full-potential  analysis  CASPOF  [Cas83]  is  used  to  provide  steady  background  flow 
information  needed  for  the  LINFLO  calculations. 

6.1  Computational  Data,  Steady  Background  Flows 

Nonlinear  steady  and  unsteady  solutions  have  been  determined  for  subsonic  and  transonic 
flows,  by  applying  the  NPHASE  analysis  on  H-type  meshes.  An  H-mesh  consisting  of  141 
axial  and  41  tangential  lines  is  used  for  the  inviscid  subsonic  calculations.  For  inviscid 
transonic  flows,  the  161  x  41  H-mesh,  shown  in  Figure  4,  is  applied,  and  141  x  81  and 
161  X  81  meshes  are  used  for  the  viscous  subsonic  and  transonic  calculations,  respectively. 
These  meshes  extend  5  axial  chords  upstream  and  9  axial  chords  downstream  from  the  blade 
row.  For  the  inviscid  subsonic  calculations,  the  average  normal  grid  spacing  adjacent  to  a 
blade  is  0.1%  of  chord  and  55  points  are  placed  along  the  blade  surfaces.  For  the  inviscid 
transonic  calculations  the  average  normal  grid  spacing  adjacent  to  a  blade  is  0.05%  of  chord, 
and  75  points  are  placed  on  the  blades.  The  mesh  spacings  for  the  viscous  calculations  are 
defined  such  that  the  law-of-the-wall  coordinate,  ,  has  an  average  value  of  approximately 
2.0  over  a  blade  surface.  This  requires  decreasing  the  average  normal  grid  spacing  adjacent 
to  the  blade  surface  to  0.005%  of  chord  for  the  viscous  subsonic  and  transonic  calculations. 
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The  full  potential  steady  (CASPOF)  and  the  linearized  inviscid  unsteady  (LINFLO) 
solutions  are  determined  on  composite  meshes  consisting  of  local  C-meshes  embedded  in 
global  H-meshes  [Cas83,  UV91].  The  C-meshes  are  used  to  accurately  resolve  the  steady  and 
unsteady  flows  around  blade  leading  edges  and  through  normal  shocks,  and  to  permit  shock 
fitting  in  the  LINFLO  calculation.  The  H-  and  C-meshes  used  in  the  LINFLO  calculations 
consist  of  155  axial  and  41  tangential  lines  and  101  radial  and  11  circumferential  lines, 
respectively.  Radial  mesh  lines  are  concentrated  near  shocks  in  the  transonic  calculations. 
A  typical  transonic  mesh  is  shown  in  Figure  5.  Since  analytic  far-field  conditions  are  applied 
in  LINFLO,  the  H-mesh  extends  only  one  axial  chord  upstream  and  downstream  from  the 
blade  row. 

The  numerical  solutions,  reported  herein,  were  determined  on  an  IBM-370  Workstation. 
Converged,  NPH.A.SE,  steady  solutions  were  obtained  after  1,000  iterations  for  an  inviscid 
calculation  and  2,000  iterations  for  a  viscous  calculation.  The  steady  subsonic  inviscid  and 
viscid  solutions  required  21  and  52  CPU  minutes,  respectively,  to  converge;  the  corresponding 
steady  transonic  solutions,  25  and  61  CPU  minutes.  The  NPHASE  unsteady  calculations 
were  started  from  the  appropriate  steady  solution,  and  performed  with  1,000  time-steps  per 
cycle  of  blade  motion  with  3  Newton  iterations  per  time  step  and  3  symmetric  Gauss-Siedel 
iterations  per  Newton  iteration.  Four  cycles  of  motion  were  used  to  converge  the  nonlinear 
inviscid  and  viscid  solutions  to  a  periodic  state.  The  subsonic  inviscid  calculations  required 
about  48  CPU  min  per  blade  passage,  and  the  subsonic  viscous  calculations  required  128 
CPU  min  per  blade  passage.  For  the  transonic  inviscid  and  viscid  solutions,  57  and  153  CPU 
min  per  blade  passage,  respectively,  were  required  for  convergence.  The  number  of  blade 
passages  required  for  a  nonlinear  unsteady  calculation  depends  upon  the  interblade  phase 
angle.  For  example,  if  a  =  60  deg,  six  passages  are  needed.  Steady  full  potential  solutions 
were  obtained,  using  CASPOF,  within  10  CPU  seconds,  and  LINFLO,  linearized,  inviscid 
solutions  required  about  90  CPU  seconds  per  unsteady  case. 

Steady  Background  Flows 

Predicted  steady  surface  pressure,  P,  and  sk-?^i-friction,  fu,,  distributions  for  the  steady 
subsonic  (M_oo  =  0.7, 17_oo  =  55  deg)  and  transonic  (M_oo  =  0.8,  fl_oo  =  58  deg)  flows 
are  shown  in  Figures  6  and  7,  respectively.  The  inviscid  (full-potential  and  Euler)  sub¬ 
sonic  predictions  for  the  blade  surface  pressures  (Figure  6a)  are  in  very  close  agreement. 
These  results  along  with  the  viscous  (thin- layer  Navier-Stokes)  predictions  indicate  that 
viscous-displacement  effects  cause  a  slight  reduction  in  the  steady  blade  pressure  loading, 
i.e.,  P(x,  0~)  —  P(x,  O"'"),  where  the  superscripts  —  and  -f-  refer  to  the  lower  (pressure)  and 
upper  (suction)  surfaces  of  the  blade,  over  the  rear  half  of  the  blade.  The  surface  skin-friction 
distributions  (Figure  6b)  indicate  that  the  viscous  flow  at  Re  =  10®  separates  (i.e.,  <  0) 

from  the  suction  surface  of  each  blade  and  reattaches  just  upstream  of  the  trailing  edge.  The 
exit  conditions,  i.e.,  the  exit  freestream  (isentropic)  Mach  number,  M+oo,  and  flow  angle, 
fl+oo,  and  the  mean  lift  force  acting  on  each  blade,  Fy  =  F  ■  e^,  as  predicted  by  the  three 
analyses  are  as  follows:  the  full-potential  indicates  that  M+oo  =  0.446,  fl+oo  =  40.2  deg,  and 
Fy  =  0.348;  the  Euler  analysis,  that  M+oo  =  0.446,  fl+oo  =  39.8  deg,  and  Fy  —  0.348;  and 
the  Navier-Stokes  analysis  yields  M+oo  =  0.461,  O+oo  =  41.2  deg  and  Fy  =  0.324.  According 
to  these  results,  viscous-displacement  effects  cause  a  7%  reduction  in  the  steady  lift  force 
that  acts  on  each  blade. 
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For  the  subsonic  cascade  configuration  operating  at  M-oo  =  0.7  and  Q-oo  =  55  deg,  if 
we  set  M+oo  =  0.446  and  =  40.2  deg,  the  resonant  interblade  phase  angles  (in  degrees) 
for  harmonic  unsteady  excitations  at  constant  interblade  phase  angle  are  =  —26.93a;, 
(7+^  =  117.12a;,  =  —31.80a;  and  =  59.79a;.  Here  the  subscripts  refer  to  the 

far  upstream  (— oo)  or  far  downstream  (+oo)  regions,  and  the  superscripts  indicate  that 
there  are  two  resonant  interblade  phase  angles  associated  with  each  of  these  regions.  The 
blade  motions  at  interblade  phase  angles  lying  between  the  lowest  (rr+oo)  and  highest 
resonant  interblade  phase  angles  are  superresonant. 

The  inviscid  surface  pressure  distributions  (Figure  7a)  for  the  steady  transonic  flow  at 
M-oo  =  0.8  and  O_oo  =  58  deg  are  in  good  agreement ,  but  differences  occur  in  the  vicinity  of 
the  shock  and  in  the  supersonic  region  lying  adjacent  to  a  blade  suction  surface  and  upstream 
of  the  shock.  Viscous  effects  reduce  the  pressure  loading  over  most  of  the  blade  and  smear 
the  pressure  change  across  the  shock.  The  transonic  flow  separates  downstream  of  the  shock 
and  is  more  or  less  separated  from  the  suction  surface  from  just  aft  of  the  shock  to  a  point 
just  upstream  of  the  trailing  edge.  For  the  transonic  flow  at  M-oo  =  0.8,  Cl-oo  =  58  deg 
the  potential  analysis  indicates  that  M^oo  =  0.432,  fl+oo  =  40.3  deg,  Fy  =  0.412,  and  the 
isentropic  Mach  numbers  at  the  base  of  each  shock  are  1.292  on  the  upstream  side  and  0.794 
on  the  downstream  side.  The  Euler  analysis  yields  M+oo  =  0.432,  fi+co  =  40.0  deg,  Fy  = 
0.408,  and  Mach  numbers  at  the  base  of  a  shock  of  1.22  and  0.871.  The  Navier-Stokes  results 
are  M+oo  =  0.461,  f2+oo  =  42.6  deg  and  Fy  =  0.372,  and  the  Mach  number  is  continuous, 
but  rapidly  varying  through  the  shock.  The  lift,  Fy,  predictions  indicate  that  viscous  effects 
lead  to  a  reduction  of  approximately  9%  in  the  lift  force  acting  on  each  blade.  For  the 
transonic  configuration  operating  at  M-oo  =  0.8  and  O_oo  =  58  deg,  if  M+oo  =  0.432  and 
0,-00  —  40.3  deg,  the  resonant  interblade  phase  angles  are;  aZoo  —  —28.94a;,  <7+^  =  201.7a;, 
(Tjioo  =  -35.92a;  and  =  66.22a;. 

6.2  Unsteady  Subsonic  Flows 

We  proceed  to  consider  unsteady  flows  through  the  subsonic  {M-oo  =  0.7,  fl_oo  =  55  deg) 
10th  Standard  Configuration  that  are  excited  by  a  unit  frequency  (a;  =  1)  in-phase  (cr  = 
0  deg)  torsional  blade  motion  with  amplitude  \a\  =  2  deg,  and  by  a  unit  frequency,  out- 
of  phase  {cr  =  180  deg),  bending  vibration  with  \hy\  =  0.01.  In  particular,  we  set  a  = 
ae^  =  (2  deg,  0)02  and  h  =  hyCy  =  (0.01,0)03,  to  define  the  torsional  and  bending  vibrations 
[cf.  (5.3)].  The  in-  and  out-of-phase  blade  motions  are  superresonant  (1,1)  and  subresonant, 
respectively.  The  instantaneous  surface  pressure  distributions  for  u>t  =90,  180,  270  and 
360  deg,  as  predicted  using  the  nonlinear  inviscid  and  viscid  {Re  =  10®)  analyses,  are  shown 
in  Figures  8  and  9,  respectively.  The  inviscid  and  viscid  unsteady  pressure  predictions  are 
very  similar.  They  indicate  that  viscous  effects  cause  a  small  reduction  in  the  time-averaged 
or  mean  pressure  loading  over  the  latter  half  of  a  blade,  but  these  effects  have  a  negligible 
impact  on  the  unsteady  surface  pressure  fluctuations. 

We  have  post-processed  the  NPHASE  subsonic,  inviscid  and  viscid,  unsteady  solutions 
to  provide  information  on  the  local  {uc)  and  global  {Wc)  works  per  cycle  in  Figures  10  and 
13,  and  on  the  harmonic  components,  p„,  of  the  unsteady  surface  pressures  in  Figures  11, 
12,  14  and  15.  The  unsteady  response  information  determined  from  the  inviscid  and  viscid 
solutions  is  nearly  identical.  Therefore,  we  include  only  the  results  of  the  nonlinear  inviscid 
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calculations  in  Figures  10  through  15,  along  with  the  linearized  inviscid  results,  determined 
using  LINFLO. 

The  predicted  local  work-per-cycle  or  pressure-displacement  function  distributions  and 
the  corresponding  global  works  per  cycle  for  the  in-phase  torsional  vibration  are  given  in 
Figure  10.  The  NPHASE  results  in  Figure  10a  are  based  on  the  inviscid  or  Euler  predictions 
for  the  nonlinear  unsteady  pressure;  those  in  Figure  10b,  on  the  Euler  predictions  for  the 
first-harmonic  component  of  the  unsteady  pressure.  The  agreement  between  the  NPHASE 
nonlinear,  the  NPHASE  first-harmonic,  and  the  LINFLO  predictions  for  and  Wc  is 
excellent,  indicating  that  the  unsteady  fluctuations  in  this  subsonic  flow  are  essentially  linear 
perturbations  of  the  underlying  steady  background  flow.  The  Navier  Stokes  solutions  for  the 
unsteady  viscous  flow  at  Rt  -  10®  provide  local  work  per  cycle  distributions  that  are  quite 
similar  to  those  given  in  Figure  10,  but  the  viscous  nonlinear  and  first-harmonic  predictions 
for  \a~'^\Wc  are  -0.61  and  -0.60,  respectively.  Thus,  the  viscous-displacement  effect  on  the 
unsteady  local  and  global  work  per  cycle  responses  is  small,  and  it  is  captured  essentially  by 
the  first-harmonic  component  of  the  unsteady  viscous  solution. 

The  real  and  imaginary  components  of  the  first  two  harmonics  of  the  unsteady  pressure 
are  shown  in  Figures  11  and  12,  respectively.  These  components  describe  unsteady  pressures 
that  are  in-(real)  and  out-of- (imaginary)  phase  with  the  blade  displacement.  The  NPHASE 
first-harmonic  inviscid  results,  shown  in  Figures  11  and  12,  as  well  as  the  corresponding  vis¬ 
cous  results  (not  shown),  are  in  very  good  agreement  with  the  LINFLO  predictions.  However, 
differences  between  the  NPHASE  and  LINFLO  predictions  for  Re{pi]  occur  on  both  the 
suction  and  pressure  surfaces  near  the  blade  leading  edge.  These  differences  could  be  due 
to  the  use  of  a  stretched  mesh  in  the  nonlinear  calculation  or,  possibly,  to  a  breakdown  in 
the  linearized  inviscid  approximation  at  blunt  blade  leading  edges.  Such  differences  do  not 
occur  in  the  results  for  Jm{|Q!|“Vi}?  nor  therefore,  in  the  local  work  per  cycle  predictions. 
The  second  (and  higher)  harmonics  of  both  the  inviscid  and  viscid  unsteady  surface  pres¬ 
sure  response  are  quite  small,  and  therefore,  these  have  a  negligible  impact  on  the  energy 
transferred  from  the  airstream  to  the  blade  {Wc)  over  a  cycle  of  blade  motion. 

Similar  results  for  the  out-of-phase  bending  vibration  are  shown  in  Figures  13  through 
15.  Again  the  NPHASE  nonlinear  and  first-harmonic  inviscid  predictions  for  the  local  and 
global  works  per  cycle  are  in  very  good  agreement  with  the  LINFLO  predictions  (Figure  13). 
The  corresponding  NPHASE  viscous  predictions  for  |Q“^|a;c  are  closely  aligned  with  the 
curves  in  Figure  13,  but,  for  the  viscous  flow  at  Re  =  10®,  the  NPHASE  nonlinear  and  first- 
harmonic  predictions  for  \a~‘^\Wc  are  both  equal  to  -8.06.  The  inviscid  unsteady  pressure 
predictions  in  Figures  14  and  15,  as  well  as  the  corresponding  viscous  predictions,  indicate 
an  excellent  agreement  between  the  NPHASE,  first-harmonic  and  the  LINFLO  predictions 
for  the  unsteady  pressure,  except  near  a  blade  leading  edge,  where  the  two  predictions  for 
the  in-phase  component  of  the  first-harmonic  unsteady  pressure  differ  significantly.  Also,  the 
second  and  higher  harmonics  of  the  unsteady  pressure  are  negligible.  Thus,  for  this  essentially 
attached,  unsteady  subsonic  flow,  the  effects  of  viscous  displacement  on  the  unsteady  pressure 
response  are  small,  and  this  response  is  essentially  captured  by  the  first-harmonic  component 
of  the  nonlinear  unsteady  pressure. 

The  behavior  of  the  global  work-per-cycle  vs.  interblade  phase  angle  is  indicated  in  Fig¬ 
ure  16  for  pure  bending  and  pure  torsional  vibrations,  at  unit  frequency,  of  the  subsonic 
10th  Standard  Cascade.  The  vertical  lines  at  the  top  of  the  figure  indicate  the  resonant 
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interblade  phase  angles  for  a  unit  frequency  excitation.  The  curves  in  Figure  16  indicate 
that  SDOF  torsional  and  bending  motions  at  unit  frequency  are  stable  {Wc  <  0)  and  that 
the  torsional  vibrations  generally  have  the  lower  stability  margin.  The  NPHASE  results 
in  Figure  16  are  based  upon  the  nonlinear  viscid  and  inviscid  solutions  for  the  unsteady 
pressure.  The  corresponding  predictions,  based  on  the  first-harmonic  component  of  the  un¬ 
steady  pressure,  are  almost  identical  to  the  results  given  in  Figure  16.  Also,  the  NPHASE 
and  LINFLO  global  response  predictions  are  in  very  good  agreement.  The  greatest  deviation 
occurs  at  cr  =  120  deg,  which  is  near  the  upstream  resonance  condition,  <j  =  The  result 
agreement  indicated  in  Figure  16  demonstrates  the  accuracy  of  both  analyses,  NPHASE  and 
LINFLO,  over  the  entire  range  of  interblade  phase  angles.  The  relatively  small  differences 
between  the  LINFLO  and  the  NPHASE  inviscid  and  viscid  solutions  indicate  that  nonlinear 
and  viscous  effects  have  a  negligible  impact  on  the  aeroelastic  stability  margins  for  pure 
torsional  and  bending  vibrations  of  the  10th  Standard  Cascade. 

The  foregoing  parametric  study  on  high  Mach  number,  subsonic,  essentially  attached, 
unsteady  flows,  excited  by  small-amplitude  blade  motions,  reveals  an  excellent  agreement 
between  the  results  of  a  nonlinear  Euler,  a  thin- layer  Navier-Stokes  and  a  linearized  inviscid 
analysis.  For  the  flows  considered,  the  essential  response  information  needed  for  aeroelastic 
studies,  is  contained  in  the  first-harmonic  component  of  the  unsteady  surface  pressure  re¬ 
sponse,  and  this  component  can  be  determined  using  a  linearized  inviscid  analysis.  Viscous- 
displacement  effects  have  very  little  impact  on  unsteady  surface  pressure  responses,  and 
therefore,  on  the  aeroelastic  stability  margin  for  a  given  blade  motion. 

6.3  Unsteady  Transonic  Flows 

Next,  we  consider  unsteady  transonic  (M_oo  =  0.8,  =  58  deg)  flows  through  the 

10th  Standard  Cascade.  Again  we  consider  blade  motions  at  a;  =  1  and  examine,  in  some 
detail,  the  unsteady  responses  to  in-phase  torsional  (at  |a|  =  1  deg)  and  out-of-phase  bending 
{\hy\  =  0.01)  vibrations.  The  torsional  amplitude  is  set  at  1  deg  so  that  the  shocks  remain 
present  throughout  the  unsteady  motion.  Our  calculations  have  shown  that  for  la]  =  2  deg, 
the  shocks  disappear  when  the  blades  are  in  their  mean  positions  (d  =  0)  and  rotating 
clockwise  at  d  =  — wa,  i.e.,  at  ujt  =  (90  -)-  360n)deg,  n  -f  0,  il,  ±2,  ....  For  the  transonic 
operating  condition  the  in-phase  blade  motions  are  superresonant  (1,1)  and  the  out-of-phase 
motions  are  superresonant  (1,0). 

Instantaneous  surface  pressure  distributions  at  tot  =  90,  180,  270  and  360  deg  as  predicted 
by  the  inviscid  and  viscous  (for  Re  =  10^)  versions  of  the  NPHASE  analysis  are  shown  in 
Figures  17  and  18,  respectively.  These  results  indicate  that  the  unsteady  pressure  fluctu¬ 
ations  on  the  suction  surface  are  considerably  stronger  than  those  on  the  pressure  surface, 
particularly  in  the  vicinity  of  the  shock.  Also,  an  observer  stationed  between  the  extreme 
shock  positions  would  experience  large  changes  in  pressure  as  the  shock  passes  by.  Finally, 
viscous  effects  tend  to  smear  the  shock,  reduce  unsteady  pressure  fluctuations  in  the  vicinity 
of  the  shock,  and  reduce  the  impulsive  loads  associated  with  shock  displacement. 

The  pressure-displacement  function  distributions  and  global  work  per  cycle  predictions 
determined  using  the  NPHASE,  thin-layer  Navier-Stokes  and  Euler  analyses  and  the  LINFLO 
linear  inviscid  analysis  are  shown  in  Figure  19  for  the  in-phase  torsional  blade  motion. 
The  NPHASE  results  in  Figure  19a  are  based  on  the  time-dependent,  nonlinear,  unsteady 
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pressure;  those  in  Figure  19b,  on  the  first-harmonic  component  of  the  unsteady  pressure. 
Shock  phenomena  are  captured  in  NPHASE  and,  although  shocks  are  resolved  accurately 
and  sharply  at  each  instant  of  time,  the  integrations  used  to  determine  the  harmonic  content 
of  the  unsteady  pressure  [cf.  (5.8)]  and  the  local  work  per  cycle  [cf.  (5.13)]  tend  to  smear 
the  unsteady  loads  in  the  vicinity  of  the  shock.  In  contrast,  LINFLO  is  a  frequency  domain 
analysis  in  which  shocks  are  “fitted”  by  imposing  shock-jump  conditions  at  the  mean  or 
steady  state  shock  positions.  This  leads  to  the  prediction  of  concentrated  shock  loads, 
represented  by  the  6  functions  in  Figure  19,  that  depend  on  the  product  of  the  steady 
pressure  jump  across  the  shock  and  the  shock  displacement,  and  first-harmonic  unsteady 
surface  pressures  that  are  discontinuous  at  the  shock,  but  continuous  elsewhere. 

The  results  in  Figure  19  indicate  that  the  nonlinear  viscid  and  inviscid  solutions  are  in 
close  agreement  except  in  the  vicinity  of  the  shock.  Viscous  displacement  effects  tend  to  shift 
the  mean  shock  position  slightly  upstream  and  to  weaken  the  impulsive  shock  load.  Also, 
local  and  global  work  per  cycle  predictions,  based  on  the  first-harmonic  component  of  the 
nonlinear  solutions  for  the  pressure,  are  in  excellent  agreement  with  the  predictions  based 
on  the  full  nonlinear  pressure.  Thus,  viscous-displacement  effects  have  only  a  small  impact, 
and  the  higher  harmonics  of  the  unsteady  pressure  response  have  a  negligible  impact,  on  the 
stability  of  the  blade  motion.  The  NPHASE  and  LINFLO  results  differ  substantially  in  the 
vicinity  of  the  shock  and  in  the  supersonic  region  upstream  of  the  shock,  but  they  are  in 
close  agreement  on  the  blade  suction  surface  downstream  of  the  shock  and  over  the  entire 
blade  pressure  surface.  Differences  at  the  shock  are  to  be  expected  because  of  the  different 
methods  used  to  determine  the  unsteady  aerodynamic  response  information.  The  differences 
in  the  supersonic  region  upstream  of  the  shock  are  unexpected  and  at  this  point  not  well 
understood.  They  could  be  associated  with  the  use  of  a  stretched  grid  for  the  NPHASE 
calculations. 

The  real  and  imaginary  parts  of  the  first  two  harmonics  of  the  unsteady  surface  pressure 
for  the  in-phase  torsional  blade  motion,  as  determined  from  the  NPHASE  viscid  and  inviscid 
solutions,  are  shown  in  Figures  20  and  21.  Similar  results  of  the  NPHASE,  nonlinear,  and 
the  LINFLO,  linearized,  inviscid  analyses  are  shown  in  Figures  22  and  23.  The  nonlinear 
viscid  and  inviscid  results  are  in  close  agreement  (Figures  20  and  21),  except  in  the  vicinity 
of  the  shock.  Also,  the  second  and  higher  harmonic  components  of  the  unsteady  pressure 
are  quite  small,  except  in  the  region  traversed  by  the  shock.  The  NPHASE  predictions  for 
the  first-harmonic  pressure  differ  substantially  from  the  corresponding  linearized  inviscid 
solutions  (Figures  22  and  23)  near  the  blade  leading  edge,  in  the  vicinity  of  the  shock,  and 
in  the  supersonic  region  upstream  of  the  shock.  The  nonlinear  and  the  linearized  inviscid 
solutions  for  the  first-harmonic  pressure  response  are  in  good  agreement  in  subsonic  regions 
of  the  flow,  i.e.,  on  the  suction  surface  downstream  of  the  shock  and  over  most  of  the  pressure 
surface.  It  should  be  noted  that  the  LINFLO  predictions  shown  in  Figures  22  and  23  only 
describe  the  first- harmonic  unsteady  pressure  and  not  the  anharmonic  unsteady  pressure 
loading  caused  by  shock  motion. 

Local  and  global  work  per  cycle  predictions  for  the  out-of-phase  bending  vibration  of  the 
transonic  cascade  are  shown  in  Figure  24.  Again,  with  the  exception  of  the  shock  loads, 
the  nonlinear  viscid  (at  Re  =  10®)  and  inviscid  solutions  are  in  good  agreement.  The 
viscous  displacement  effect  shifts  the  mean  shock  location  slightly  upstream  and  weakens 
the  unsteady  shock  loads.  The  results  for  ojc{x)  and  Wc,  based  on  the  nonlinear  and  the 
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first-harmonic  unsteady  pressure,  are  in  very  close  agreement  indicating  that  the  energy 
transfer  between  the  fluid  and  the  blade  motion  is  almost  entirely  due  to  the  first-harmonic 
response  characteristics  of  the  fluid.  The  linearized  inviscid  analysis,  LINFLO,  predicts  very 
high  and  very  low  values  of  the  local  work  per  cycle  just  upstream  and  just  downstream  of 
the  shock,  respectively.  These  are  caused  by  the  first-harmonic  unsteady  pressure  response. 
The  LINFLO  analysis  also  predicts  a  laxge  anharmonic  local  unsteady  response  due  to  shock 
motion,  which  is  represented  by  the  delta  functions  in  Figure  24.  The  NPHASE  and  LINFLO 
results  on  the  suction  surface,  away  from  the  shock,  and  on  the  pressure  surface  are  in  good 
agreement,  but  the  two  analyses  predict  very  different  shock  loads. 

The  NPHASE  viscid  and  inviscid  predictions  for  the  real  and  imaginary  components  of 
the  first  and  second  harmonics  of  the  unsteady  surface  pressure,  for  the  out-of-phase  bending 
vibration  are  shown  in  Figures  25  and  26.  The  corresponding  NPHASE  nonlinear  inviscid 
and  LINFLO  linearized  inviscid  results  are  given  in  Figures  27  and  28.  Again,  the  results 
of  the  nonlinear  viscid  and  inviscid  analyses  are  in  close  agreement,  except  in  the  region 
traversed  by  the  shock,  where  a  relatively  large  second  and  higher  (not  shown)  harmonic 
components  of  the  unsteady  pressure  persist.  The  NPHASE  inviscid  predictions  for  the  real 
and  imaginary  components  of  the  first-harmonic  unsteady  pressure  differ  substantially  from 
the  LINFLO  predictions  at  the  blade  leading  edge,  in  the  supersonic  region  upstream  of  the 
shock,  and  in  the  vicinity  of  the  shock. 

Global  work-per-cycle  predictions,  based  on  the  NPHASE  Navier-Stokes  and  Euler  analy¬ 
ses  and  on  the  LINFLO  linearized  inviscid  analysis,  are  shown  in  Figure  29  for  unit-frequency 
torsional  and  bending  vibrations  of  the  10th  Standard  Cascade.  The  NPHASE  viscid  and 
inviscid  results  are  in  fairly  close  agreement,  indicating  that  for  the  unsteady  flows  considered 
viscous-displacement  effects  are  of  limited  importance,  but  they  do  tend  to  reduce  the  stabil¬ 
ity  margins  for  the  torsional  blade  vibrations.  The  results  of  the  NPHASE  analyses  and  the 
LINFLO  analysis  show  similar  trends,  but  significant  quantitative  differences,  particularly 
for  the  bending  vibrations.  These  differences  are  due  to  the  relatively  large  discrepancies  be¬ 
tween  the  local  response  predictions  of  the  two  analyses  at  shocks  and  in  supersonic  regions. 
Work  per  cycle  predictions,  based  on  the  NPHASE  viscid  and  inviscid  predictions  for  the 
first-harmonic  unsteady  pressure,  agree  very  closely  with  the  nonlinear  results  shown  in  Fig¬ 
ure  29,  again  indicating  that,  even  in  transonic  flows,  the  stability  margin  for  a  given  blade 
motion  is  essentially  determined  by  the  first-harmonic  component  of  the  unsteady  pressure 
response. 

The  results  of  the  foregoing  parametric  study  for  transonic  flows  are  not  as  comforting  or 
convenient  as  those  determined  previously  for  subsonic  flows.  The  most  important  concerns 
are  the  discrepancies  between  the  nonlinear  predictions  and  the  linearized  inviscid  predictions 
in  supersonic  regions  and  in  the  vicinity  of  shocks.  The  discrepancies  in  supersonic  regions 
could,  perhaps,  be  eliminated  by  incorporating  accurate,  unsteady,  far-field  conditions  into 
the  NPHASE  analysis,  to  avoid  the  use  of  stretched  H-meshes  (see  Figure  4).  The  ways  in 
which  shocks  and  their  motions  are  modeled  in  nonlinear  and  linearized  unsteady  codes  must 
be  understood  more  clearly  to  resolve  the  differences  between  the  NPHASE  and  LINFLO 
predictions  in  the  vicinity  of  shocks.  Finally,  although  viscous-displacement  effects  are  gen¬ 
erally  small  for  unsteady  transonic  flows  with  small-scale  boundary  layer  separations,  such 
effects  can  be  important  shock  regions.  Therefore,  viscous  effects  should  be  incorporated  into 
the  approximate  unsteady  aerodynamic  analyses  intended  for  transonic  aeroelastic  design  ap- 
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plications.  An  encouraging  result  of  the  transonic  parametric  studies  is  the  close  agreement 
obtained  between  the  local  and  global  work  per  cycle  predictions  that  were  based  on  the 
NPHASE  nonlinear  and  first-harmonic  pressure  responses.  This  suggests  that  a  linearized 
analysis  that  accounts  for  viscous-displacement  effects  would  be  useful  for  predicting  the 
flutter  characteristics  of  turbomachinery  blading  over  a  wide  range  of  operating  conditions. 
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7.  Concluding  Remarks 


A  modern  implicit,  time-marching,  Euler/Navier-Stokes  analysis  (NPHASE)  has  been  ap¬ 
plied  to  simulate  two-dimensional  subsonic  and  transonic  unsteady  flows  through  a  vibrating 
cascade,  in  particular,  the  Tenth  (10th)  Standard  Cascade  of  Ref.  [FV93].  The  vibrations 
considered  are  pure  translations  normal  to  the  blade  chords  (bending)  and  pure  rotations 
about  axes  at  blade  midchords  (torsion).  Linearized  inviscid  response  predictions,  as  deter¬ 
mined  using  the  LINFLO  code,  have  also  been  determined  to  help  validate  and  illustrate 
the  Euler/Navier-Stokes  capabilities  for  unsteady  cascade  flows.  The  NPHASE  numerical 
results  have  also  been  interrogated  to  assess  the  importance  of  nonlinear  and  viscous  effects 
on  the  unsteady  aerodynamic  response  of  blades  vibrating  at  a  frequency  representative  of 
those  at  which  flutter  occurs  in  modern  fans  and  compressors. 

The  numerical  simulations  for  the  steady  and  unsteady  subsonic  flows  indicate  that  vis¬ 
cous  effects  tend  to  have  a  small  impact  on  steady  blade  loading  and  a  negligible  impact  on 
the  unsteady  loading.  In  addition,  the  unsteady  surface  pressure  fluctuations  are  essentially 
linear.  The  results  of  the  NPHASE  Euler  analysis  and  a  full  potential  analysis  for  steady 
flow  are  in  excellent  agreement.  The  steady  viscous  flow  predictions,  for  Re  =  10®,  exhibit  a 
small  decrease  in  the  blade  pressure  loading  on  the  pressure  surface  over  the  rear  half  of  the 
blade  and  a  small  region  of  separated  flow  on  the  suction  surface  extending  from  88  to  98%  of 
the  blade  chord.  The  unsteady  simulations  show  that  the  viscid  and  inviscid  response  predic¬ 
tions,  i.e.,  the  local  and  global  works  per  cycle  and  the  unsteady  surface  pressures  are  in  close 
agreement.  The  stability  of  the  blade  motion  is  essentially  determined  by  the  first-harmonic 
component  of  the  unsteady  pressure,  as  the  second-  and  higher-harmonic  components  are 
quite  small.  The  LINFLO  linearized  inviscid  predictions  for  the  unsteady  response  quantities 
are  in  excellent  agreement  with  both  the  NPHASE  Euler  and  Navier-Stokes  predictions. 

.  For  the  steady  and  unsteady  transonic  flows,  the  numerical  simulations  indicate  that 
accurate  predictions  of  steady  and  unsteady  surface  pressures  depend  significantly  upon  an 
adequate  resolution  of  shock  strength  and  the  positions  at  which  shocks  impinge  on  blade 
surfaces.  For  the  steady  flow,  the  Euler  and  Full  Potential  results  for  the  surface  pressure 
are  in  excellent  agreement  at  points  along  the  blade  surface  at  which  the  flow  is  subsonic, 
but  differ  somewhat  in  the  supersonic  region  upstream  of  the  suction-surface  shock  and  in 
the  vicinity  of  this  shock.  Navier-Stokes  predictions  for  the  corresponding  viscous  flow  at 
i?e  =  10®  indicate  that  viscous  displacement  effects  tend  to  weaken  and  smear  the  shocks, 
increase  the  pressures  on  the  suction  surface  in  the  vicinity  of  the  leading  edge  and  reduce 
the  pressures  along  the  pressure  surface  of  the  blade.  The  viscous  steady  flow  separates 
from  the  suction  surface  downstream  of  the  shock,  and  is  either  separated  or  on  the  verge  of 
separation  from  the  shock  to  a  point  just  upstream  of  the  trailing  edge. 

Predictions  for  the  unsteady  transonic  flows  indicate  that  the  unsteady  pressure  fluctu¬ 
ations  on  the  suction  surface  are  much  stronger  than  those  on  the  pressure  surface,  with 
severe  fluctuations  in  surface  pressure  resulting  from  the  motion  of  a  shock.  Viscous  ef¬ 
fects  diminish  the  strength  of  the  impulsive  unsteady  loads  associated  with  shock  motion. 
Also,  viscous-displacement  effects  have  only  a  small  impact  on  global  and  local  work  per 
cycle  predictions,  causing  significant  changes  in  the  local  work  per  cycle  predictions  only 
in  the  vicinity  of  a  shock.  The  transonic  unsteady  simulations  indicate  that  large  second- 
and  higher-harmonic  pressure  variations  occur  in  the  region  traversed  by  a  shock,  but  the 
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higher-harmonic  content  of  the  unsteady  pressure  response  has  only  a  negligible  impact  on 
the  local  and  global  work  per  cycle  responses,  and  hence,  on  the  stability  margins  for  the 
blade  motions.  Unfortunately,  the  LINFLO  linearized  inviscid  results  are  not  in  good  agree¬ 
ment  with  the  nonlinear  predictions  in  the  vicinity  of  shock  and  in  the  supersonic  region 
upstream  of  the  shock.  The  linearized  inviscid  results  are  in  reasonable  agreement  with  the 
nonlinear  inviscid  and  viscid  predictions  in  subsonic  regions  of  the  flow. 

There  remain  several  issues  which  require  further  study  before  the  impact  of  nonlinear 
and  viscous  displacement  effects,  as  they  apply  to  blade  flutter  prediction,  can  be  fully 
understood.  The  present  investigation  has  indicated  that  fully  two-dimensional,  unsteady 
far-field  conditions  should  be  implemented  into  the  NPHASE  code  to  avoid  the  need  for 
calculating  unsteady  flows  on  stretched  H-meshes.  If  NPHASE  and  LINFLO  calculations  are 
performed  on  similar  H-meshes,  the  discrepancies  between  the  results  of  the  two  analyses  in 
supersonic  regions  and  near  blade  leading  edges  may  disappear.  Also,  a  better  understanding 
of  shock  phenomena  is  warranted.  The  goals  here  should  be  to  bring  the  NPHASE  and 
LINFLO  unsteady  transonic  predictions  into  closer  agreement,  and  to  understand  more 
clearly  the  reasons  for  the  differences  between  the  shock  loads  predicted  by  the  two  analyses. 

The  parametric  studies  conducted  under  the  present  effort  pertained  to  inviscid  flows 
and  fully-turbulent,  viscous  flows  at  high  Reynolds  number.  Any  shocks  occurring  in  these 
flows  were  present  throughout  the  unsteady  motion,  and  boundary  layer  separations  were  of 
limited  extent.  NPHASE  has  proven  to  be  a  robust  code  for  such  applications,  which  can  give 
useful  insights  into  the  effects  of  nonlinear  and  viscous  unsteady  aerodynamic  phenomena 
on  blade  response.  Future  research  should  be  directed  towards  investigating  unsteady  flows 
in  which  shocks  appear  and  disappear  within  a  cycle  of  blade  motion  to  gain  insight  into  the 
importance  of  large  shock  excursions  and  their  impact  on  aeroelastic  stability.  Also,  unsteady 
flows  in  which  large-scale  viscous-layer  separations  occur  should  be  investigated.  Finally, 
the  impact  of  various  turbulence  and  transition  models  on  predicted  unsteady  aerodynamic 
response  behaviors  should  also  be  investigated. 
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Figure  2:  Unsteady  excitations:  blade  motion,  incident  vortical  and  entropic  disturbances 
from  upstream,  and  incident  acoustic  disturbances  from  upstream  and  downstream. 
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Figure  4:  H-mesh  used  for  the  NPHASE,  steady  and  unsteady,  inviscid,  transonic  simula¬ 
tions. 
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Figure  5:  Composite  H/C  mesh  used  for  the  LINFLO  unsteady  transonic  calculations. 
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Figure  6;  Predicted  steady  surface  pressure  (a)  and  skin  friction  (b)  distributions  for  inviscid 
and  viscid  (at  Re  —  10®)  subsonic  (M_oo  =  0.7,  fi-oo  =  55  deg)  flows  through  the  10th 
Standard  Cascade. 
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_ Euler  Analysis  (NPHASE) 
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Figure  7:  Predicted  steady  surface  pressure  (a)  and  skin  friction  (b)  distributions  for  inviscid 
and  viscid  (at  Re  =  10®)  transonic  {M-^o  -  0.8,  =  58  deg)  flows  through  the  10th 

Standard  Cascade. 
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Figure  8:  Instantaneous  surface  pressure  distributions  for  inviscid  flows  through  the  subsonic 
(M_oo  =  0.7,  i7_oo  =  55  deg)  10th  Standard  Cascade  undergoing  prescribed  blade  vibrations: 
(a)  in-phase  torsional  vibration  about  midchord  at  |q:|  =  2  deg,  u)  =  1,  and  a  =  0  deg;  (b) 
out-of-phase  bending  vibration  at  l/iy]  =  0.01,u;  =  1,  a  =  180  deg. 
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Figure  9:  Instantaneous  surface  pressure  distributions  for  viscous  flows  at  Re  =  10®  through 
the  subsonic  (JV/_oo  =  0.7,  fl-oo  =  55  deg)  10th  Standard  Cascade  undergoing  prescribed 
blade  vibrations:  (a)  in-phase  torsional  vibration  about  midchord  at  |a|  =  2  deg,  a;  =  1,  and 
(7  =  0  deg;  (b)  out-of-phase  bending  vibration  at  |/iy|  =  0.01,  u  =  1,  a  =  180  deg. 
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Figure  10:  Pressure-displacement  function  distributions  for  inviscid  subsonic  flow  (M_oo  == 
0.7,  Q-oo  =  55  deg)  through  the  10th  Standard  Cascade  undergoing  an  in-phase  torsional 
vibration  about  midchord  at  |a|  =  2  deg,  u>  =  1,  and  a  =  0  deg:  (a)  NPHASE  nonlinear 
prediction;  (b)  NPHASE  first-harmonic  prediction. 
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Figure  11:  Real  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure  dis¬ 
tribution  for  inviscid  subsonic  (M_co  =  0.7,  Vl-oo  =  55  deg)  flow  through  the  10th  Standard 
Cascade  undergoing  an  in-phase  torsional  vibration  about  midchord  at  |q:|  =  2  deg,  a;  =  1, 
and  (j  =  0  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  12:  Imaginary  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  inviscid  subsonic  (M_oo  =  0.7,  fi_oo  =  55  deg)  flow  through  the  10th  Stan¬ 
dard  Cascade  undergoing  an  in-phase  torsional  vibration  about  midchord  at  [a]  =  2  deg, 
w  =  1,  and  a  =  0  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  13:  Pressure-displacement  function  distribution  for  inviscid  subsonic  (M_oo  =  0.7, 
f^-oo  =  55  deg)  flow  through  the  10th  Standard  Cascade  undergoing  an  out-of-phase  bending 
vibration  at  \hy\  =  0.01,  tu  =  1,  and  a  =  180  deg:  (a)  NPHASE  nonlinear  prediction;  (b) 
NPHASE  first-harmonic  prediction. 
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Figure  14:  Real  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  inviscid  subsonic  {M-oo  =  0.7,  fi_oo  =  55  deg)  flow  through  the  10th  Stan¬ 
dard  Cascade  undergoing  an  out-of-phase  bending  vibration  at  |hj,|  =  0.01,  a;  =  1,  and 
(7  =  180  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  15:  Imaginary  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pres¬ 
sure  distribution  for  inviscid  subsonic  —  0.7,  fi-oo  =  55  deg)  flow  through  the  10th 

Standard  Cascade  undergoing  an  out-of-phase  bending  vibration  at  \hy\  =  0.01,  w  =  1,  and 
cr  =  180  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  16:  Work  per  cycle  versus  interblade  phase  angle  for  inviscid  and  viscous  (at  Re  =  10®) 
subsonic  flows  through  the  10th  Standard  Cascade  undergoing  prescribed  blade  motions  at 
unit  frequency  (a;  =  1):  (a)torsional  vibrations  at  |q;|  =  2  deg  about  midchord;  (b)  bending 
vibrations  at  |/iy|  =  0.01. 
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Figure  17:  Instantaneous  surface  pressure  distributions  for  inviscid  flows  through  the  tran¬ 
sonic  (M_oo  =  0.8,  Cl-oo  =  58  deg)  10th  Standard  Cascade  undergoing  prescribed  blade 
vibrations:  (a)  in-phase  torsional  vibration  about  midchord  at  |q:|  =  1  deg,  w  =  1,  and 
(7  =  0  deg;  (b)  out-of-phase  bending  vibration  at  \hy  \  —  0.01,  a;  =  1,  (7  =  180  deg. 
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Figure  18:  Instantaneous  surface  pressure  distributions  for  viscous  flows  (at  Re  =  10®) 
through  the  transonic  (M_oo  =  0.8,  Cl-oo  =  58  deg)  10th  Standard  Cascade  undergoing 
prescribed  blade  vibrations:  (a)  in-phase  torsional  vibration  about  midchord  at  |q:|  =  1  deg, 
u)  =  1,  and  (7  =  0  deg;  (b)  out-of-phase  bending  vibration  at  |fty|  =  0.01,  u>  =  1,  <t  =  180  deg. 
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Figure  19:  Pressure-displacement  function  distribution  for  viscid  (at  Re  =  10®)  and  inviscid 
transonic  flows  {M-oo  =  0.8,  O_oo  =  58  deg)  through  the  10th  Standard  Cascade  undergoing 
an  in-phase  torsional  vibration  about  midchord  at  |q;|  =  1  deg,  a;  =  1,  and  (7  =  0  deg:  (a) 
NPHASE  nonlinear  predictions;  (b)  NPHASE  first-harmonic  predictions. 
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Figure  20:  Real  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  viscid  (at  Re  =  10®)  and  inviscid  transonic  (M_oo  =  0.8,  Q-oo  =  58  deg) 
flows  through  the  10th  Standard  Cascade  undergoing  an  in-phase  torsional  vibration  about 
midchord  at  |a|  =  1  deg,  u;  =  1,  and  cr  =  0  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  21:  Imaginary  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  viscid  (at  Re  =  10®)  and  inviscid  transonic  (M_oo  =  0.8,  H-oo  =  58  deg) 
flows  through  the  10th  Standard  Cascade  undergoing  an  in-phase  torsional  vibration  about 
midchord  at  |a|  =  1  deg,  a>  =  1,  and  <7  =  0  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  22:  Real  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  inviscid  transonic  {M-oo  =  O-Sj  f^-oo  =  58  deg)  flow  through  the  10th 
Standard  Cascade  undergoing  an  in-phase  torsional  vibration  about  midchord  at  |a|  =  1  deg, 
a;  =  1,  and  (T  =  0  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  23:  Imaginary  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  inviscid  transonic  {M-^  =  0.8,  =  58  deg)  flow  through  the  10th 

Standard  Cascade  undergoing  an  in-phase  torsional  vibration  about  midchord  at  |q:|  =  1  deg, 
a;  =  1,  and  a  =  0  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  24:  Pressure-displacement  function  distribution  for  viscid  (at  Re  —  10®)  and  inviscid 
transonic  (M_co  =  0.8,  fl_oo  =  58  deg)  flows  through  the  10th  Standard  Cascade  undergoing 
an  out-of-phase  bending  vibration  at  =  0.01,  w  =  1,  and  cr  =  180  deg:  (a)  NPHASE 
nonlinear  predictions;  (b)  NPHASE  first-harmonic  predictions. 
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Figure  25:  Real  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  viscid  (at  Re  =  10^)  and  inviscid  transonic  (M^oo  =  0.8,  fi_oo  =  58  deg) 
flows  through  the  10th  Standard  Cascade  undergoing  an  out-of-phase  bending  vibration  at 
|/it/|  =  0.01,  a;  =  1,  and  cr  =  180  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  26:  Imaginary  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  viscid  (at  Rt  =  10®)  and  inviscid  transonic  (M_oo  =  0.8,  f^-oo  =  58  deg) 
fiows  through  the  10th  Standard  Cascade  undergoing  an  out-of- phase  bending  vibration  at 
=  0.01,  w  =  1,  and  a  =  180  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  27:  Real  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  inviscid  transonic  (Af_oo  =  0.8,  O_oo  =  58  deg)  flow  through  the  10th 
Standard  Cascade  undergoing  an  out-of-phase  bending  vibration  at  \hy\  =  0.01,  a;  =  1,  and 
cr  =  180  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  28:  Imaginary  components  of  the  Fourier  amplitudes  of  the  unsteady  surface  pressure 
distribution  for  inviscid  (M_oo  =  0.8,  fl-oo  =  58  deg)  flow  through  the  10th  Standard 
Cascade  undergoing  an  out-of-phase  bending  vibration  at  =  0.01,  u;  =  1,  and  <7  = 
180  deg:  (a)  Suction  surface;  (b)  Pressure  surface. 
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Figure  29:  Work  per  cycle  versus  interblade  phase  angle  for  inviscid  and  viscous  (at  Re  =  10®) 
transonic  flows  through  the  10th  Standard  Cascade  undergoing  prescribed  blade  motions  at 
unit  frequency  (a;  =  1):  (a)  torsional  vibrations  at  |q;|  =  1  deg  about  midchord;  (b)  bending 
vibrations  at  |Ay|  =  0.01. 
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